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ABSTRACT 


A theory of radial orbit errors is derived from first principles; the 
basic ideas are presented in a reasonably self-contained way. There is 
also a section about tides, their effects on orbits, and their possible 
mapping with satellite altimetry. The theory seems to explain well the 
outcomes of computer simulations made independently by this author and by 
others. Orbits with precisely repeating ground tracks are considered in 
detail. SEASAT was kept in an orbit like that during its last month of 
operation, and future satellites with altimeters are likely to be in 
orbits of this type (TOPEX, ERS-1, POSEIDON). In this special case, 

Ar^®\ the part of the radial error caused by the gravitational field 

model, has a component that depends only on position along the repeating 
ground track; the remainder of A r is caused by deep orbital resonance 
(mostly with the zonals) and consists mainly of a slowly increasing 
oscillation of one cycle per revolution (non-gravitational force models 
may produce similar "resonant" errors). Moreover, the periodical errors 
due to the zonals are functions of latitude only, and unobservable in 
crossover differences. Repeat orbits, where the perigee librates very 
slowly, are not the most suitable, therefore, for using altimetry to 

compute the mean sea surface by the "bi&s+tilt" and similar methods, or 

for gaining new information about the zonal part of the field. These 
orbits are best for studying temporal changes of the sea surface, because 
the error can be filtered out easily using colinear pass differences. 
Conversely, non-repeating orbits, where the perigee precesses, are worse 
for studying changes but are better for mapping the mean sea surface and 
modelling the field, provided the altimeter observations span a 
substantial part of one apsidal cycle. 
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1.0 INTRODUCTION 

The face of the ocean has been examined with artificial satellites in 
the last decade and found to be covered with remarkable features hardly 
suspected before, when more than two thirds of the globe were still 
largely unsurveyed. Among the millions of measurements made with diverse 
instruments, those taken with the radar altimeters of GEOS 3 and SEASAT 
have proven of particular value for studying our planet, because of their 
high and homogeneous quality, dense coverage, and range of application. 
So it is not surprising that new space missions involving the use of 
altimeters are being planned for the next decade: those of TOPEX (NASA), 

ERS-1 (ESA) and POSEIDON (ONES), among others. 

The importance- of altimetry is suggested by the volume and variety of 
the scientific literature related to it. One could mention, for example, 
three special Issues of the Journal of Geophysical Research (Vol. 84, B8, 
1979 on GEOS 3; and Vols. 87, C5, 1982 and 88, C3, 1983 on SEASAT), one of 
Marine Geodesy (Vol. 8, 1984), or one of the Journal of the Astronautical 
Sciences (Vol. 28, October-December 1980 on orbit determination for 
SEASAT), as well as the proceedings of a number of international 
symposia. Two comprehensive reviews of the work published in the USA 
alone from 1975 to 1982 have been made by Stanley ( 1979) and by Marsh 
(1983). 
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1.1 Motivation 


Since the earliest days of GEOS 3 it has been clear that the orbit 
error constitutes a serious problem. Things have improved, but the 
problem persists. While the error has its uses, too, as a source of 
information on what is not known about the gravitational field (Shum, 
1982, and Wagner, 1984), it remains mostly a nuisance that can be 

mitigated with more or less empirical remedies, but that still limits the 
usefulness of the data in many applications. For this reason, efforts to 
get better models of the forces acting on the spacecraft, of which the 

gravitational ones arfi the most significant, are currently under way on 
both 9ides of the Atlantic. 

An approach to the orbit problem based on celestial mechanics can 
give further insight into the way the empirical methods work and, perhaps, 
into how to improve them. As the relevant literature on satellite geodesy 
includes many important papers, reports and books either out of print or 
hard to get, I have written here a rather detailed introduction to the 
basic concepts. This may be of some value to those who wish to understand 
the nature of the orbit error better and are not very familiar with the 
subject. 

The theory, as given here, is slanted towards orbits determined by 

the methods now in common use for geodetic satellites (Martin et al. , 

1976; Putney, 1976; Schutz et al., 1980). The orbits of GEOS 3 and SEASAT 
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have been compiled by those methods, which are also the main ones 
considered, at this stage, for future missions. One possible alternative 
is to find the position of the spacecraft from simultaneous observations 
of several GPS-NAVSTAR navigation satellites. If the ephemerldes of these 
satellites were sufficiently well known, this would eliminate the need to 
integrate long arcs numerically with imperfect force models, thus avoiding 
the orbit errors due to those models. The idea has been studied by Bender 
and Larden (1982) and by Ondrasik and Wu (1982). It may be used, on an 
experimental basis, during part of the TOPEX mission. 

The analytical orbit perturbation theory explained here is only meant 
to give a qualitative understanding of the error. Its formulation 
Includes several simplifying assumptions, so it is quantitatively less 
accurate than the numerical integration of the exact differential 
equations Involved, which is also much more adequate for handling 
non-gravitational forces (drag, radiation pressure, etc.), and so it has 
become the standard procedure in satellite geodesy today (its main 
limitation complementing that of the analytical approach: it trades 

insight for numbers). Nevertheless, efforts to refine and extend the 
analytical method still continue, as shown in a review by Gaposchkin 
(1978). The treatment of the radial error along a precisely repeating 
ground track is a new version, using Keplerian elements, of part of a 
complete theory developed by me for this special type of orbit while 
visiting at the Department of Geodesy of the Technical University at 
Delft, in the Netherlands, over the period 1982-83 (Colombo, 1984). More 
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recently, a similar theory for the periodical part of that error has been 
worked out independently by C.A. Wagner (ib., 1984). 

While oriented towards oceanography (the main field of application of 
altimetry until now) the present discussion of the orbit error is also 
relevant to the uses of data taken on land, ice, etc. 

The topic of tides is intimately linked to that of altimetry. Tides 
modify the gravitational field, affecting the orbits of spacecraft and 
contributing to the errors in their determination. They also show up as 

part of the signal in the altimeter measurements. In principle, it is 
possible to map them directly over the whole of the oceans using 

tjLtimetry. Attempts have been made several times, mostly with scant 

success. This must be blamed, to some extent, on the inadequate temporal 
and spatial sampling and limited accuracy of the GEOS 3 data, and on the 
short life of SEASAT. But, perhaps, better results could be obtained, 
even with such measurements, if one were to take full advantage of what 
can be predicted from first principles about the nature of the orbit 
errors. As mentioned in the section on "Tides", some of the best results 
so far have come from data collected during the last useful month of 

SEASAT. At that time, the satellite had a precisely repeating ground 

track which, in theory, should make the radial error particularly easy to 
filter, because of the strong symmetry of the orbit. 
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The literature on tides is very extensive. Cartwright estimated in 
1977 that, over the previous 200-odd years, the number of publications had 
grown almost exponentially with time (to Illustrate this , he tabulated the 
increment in this number at fifty years' intervals). Hendershott and Munk 
(1970), Cartwright (1977) and Schwiderski (1980) have made comprehensive 
reviews. The basic ideas needed to consider the tides in the context of 
satellite altimetry are scattered far and wide in works on astronomy, 
geodesy, oceanography and geophysics. I have tried to bring these ideas 
together in a way that shows, as clearly as possible, the links between 
the theory of tides and that of orbits. 

Tides have always fascinated those with imagination: regular 
movements in an otherwise unpredictable element, following the rhythms of 
the sky. Now they are beginning to be mapped from that same sky, with 
modern instruments whose proper use requires understanding the timeless 
rules by which the sea, the Sun, the Moon, the satellites and the Earth 
join together in a complex dance. In the end, through the combination of 
techniques old and new, a task formulated nearly three centuries ago shall 
be fulfilled. At the very dawn of modern celestial mechanics and tidal 
theory, Newton wrote in his "Principia": 

"Thus have I explained the causes of the motions of the 
Moon and the sea. Now it is fit to subjoin something 
concerning the amount of those motions." 
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2.0 ALTIMETRY 


2. 1 Residual Sea Heights 


The height h w of the sea surface above a reference ellipsoid , or 
ellipsoidal height of the surface, can be expressed as the sum 


h w " N + w o + C G + w t 


( 1 ) 


where N is the geoid height, w Q is the constant part of the difference 
between hy and N, or stationary sea surface topography , represents the 
total geocentric tide (a regular movement of the sea surface towards and 
away from the center of mass of the Earth, or geocenter) , while w t 
corresponds to all temporal variations in surface height other than tides. 


The value of N can be calculated using a gravitational field model. 
This model, for some ocean studies and for the computation of satellite 
orbits, has usually the form of a sum of spherical harmonics. In this 
work I shall assume that the field model is of this type. The high 
frequency part of the geoid at sea is known quite well nowadays thanks to 
the two altinetric missions already carried out, GEOS-3 and SEASAT. 
Spherical harmonic models up to degree and order 180 (and even higher) 
have been obtained from combinations of altimetry with land gravimetry, 
and they appear to be quite reliable over the oceans, particularly in 
details smaller than 1000 km (see, for example, Lerch et al. , 1981, and 
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Rapp, 1982). As shown by calibration against independent data and each 
other, existing models seem to be weakest in the frequency range between 
degrees 10 and 40. Information on the low degrees comes mostly from 
satellite tracking data. Recent intercomparisons of satellite-derived 
nodels have been made by Reigber (1983), Lambeck and Coleman (1983), and 
Lerch et al. (1984), The broadest features (below degree 8) may be quite 
reliable already (Wagner, 1983). 

The contemporary ocean tidal charts are believed to be reasonably 
good (Schwiderski, ib. 1980), but direct evidence for this is limited to 
data from tidal stations scattered widely along the coasts and in islands, 
most of which have been used to make the charts in the first place. Solid 
earth tides are thought to be known much better than ocean tides at 
present. No model for the prediction of w t exists today, though some of 
the phenomena involved are understood to some extent, such as raesoscale 
eddies, the piling up of water along the coast caused by wind, etc. There 
is reason to believe that much of w t can be treated, geographically, 
mostly as a random variable (see Wunsch, 1980, for further details) whose 
standard deviation, while changing from place to place according to how 
energetic are the currents and winds that agitate the sea surface, can be 
estimated directly from the study of overlapping altimetry passes as in 
(Cheney et al. , 1983), for example. Maps of the mean sea surface 
topography, based on general circulation models, are also available, but 
their accuracy is unclear. As examples one can mention those by Lisitzin 
(1974), Levitus and Dort (1977), etc. Attempts to map w Q using altimetry 
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have been severely hampered by the orbit and geoid errors, directly linked 
to the long wave uncertainties in existing gravitational field models. 
This problem was recognized very long ago, as shown in a paper written by 
Von Arx (1966) almost a decade before the launch of GEOS 3. Wunsch and 
Gaposchkin (1980) have proposed the simultaneous determination of w Q and N 
by analyzing a combination of gravitational, altimetric and oceanographic 
data with a suitable form of generalized linear regression estimator (see 
Rao, 1965, Ch. 2, (g) and also Moritz (1980)). Extensive discussions of 
the sea surface topography problem, and reviews of much of the work done 
using GEOS 3 data, have been written by Rizos (1980) and Coleman (1981). 

If N c is the geoid height computed from the field model, and AN the 
corresponding error in N c , then 

N - N - AN . (2a) 

c 

Similarly, if A? G is the error in ? Gc> the geocentric tide calculated 
from charts of the ocean and the Earth tide with the corresponding Love 
numbers (expression (77) in the section "Tides"), 

<G' Sjc ; ( 2 W; 
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If h_ is the ellipsoidal height of the altimeter satellite, and A the 
s 

shortest distance between the center of mass of the spacecraft and the sea 
surface, then, as a very good approximation, 

h » h - A . (3) 

w s 

After several corrections have been made to the raw data ( "pre- 
processing”, see Hancock et al. , 1980, and Tapley et al. , 1982), including 
those resulting from "in-flight" calibration of the altimeter 
(Kolenkiewicz and Martin, 1982; Marsh and Williamson, 1982), the radar 
altimeter ranges A^ can be regarded as successive values of A plus a 
measurement error AA, where AA is mostly random noise, For GEOS-3, the 
first satellite dedicated to altimetric observations , AA stood at about 
75 cm; the error came down to near 10 cm for SEASAT, launched three 
years later (in 1978); for future missions, a third generation of 
altimeters is likely to bring AA Into the centimeters' range (see 
MacArthur, 1980) » 

The height of the satellite h Q can be calculated from its computed 
orbit , or ephemeris, which is obtained by adjusting the initial position, 
velocity and a few other parameters in an iterative Gauss-Newton (or 
similar) procedure to minimize the . sum of the squares of the differences 
between the values of the tracking data available and the corresponding 
values calculated from the estimated ephemeris (see Martin et al., ib., 
1976). This fitted orbit usually spans a few days (up to one week), and 
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then a new orbit fit is made for the following several days, and so on. 

While not unique, this is today a common procedure for getting accurate 

ephemerldes of altimeter satellites, so in this study I consider only 

orbits computed in this way. At present, the error Ah in h is believed 

s s 

to be of the order of 1.5 meters (r.m.s.) (see Marsh and Williamson, 1980; 
Lerch et al. , 1982). CD 

The reasons Why the ephemsris are not exact are multiple: errors in 

tracking data and in tracking station coordinates; the abundance and 

distribution, temporal as well as geographical, of these data; imperfect 
refraction corrections in the tracking; errors in the models for 

calculating the forces that shape the orbit, etc. Of the force errors, 
the most prominent, at present, are errors in the gravitational field 
model, followed at some distance by errors in the models of surface forces 
such as air drag and solar radiation pressure. The last two are 

particularly serious in satellites like SEASAT, because of the large solar 
panels needed to feed their power-hungry equipment, as well as the 

presence of a number of large antennas and other objects of complex shape, 
all of which is made worse by the changing attitude of the spacecraft, as 
this varies the effective cross-section opposed to the surface forces. 
Often today orbits are estimated assuming that u satellite is a 

^ ^Horizontal orbital errors can be disregarded. They are, at present, 

of the order of less than 100 m, and because of the smoothness of the sea- 
surface (averaged over the crost-section of the radar altimeter beam, or 
footprint, usually 1 km or more in diameter) their effects on the accuracy 
of the calculated height is negligible. Numerical integration errors in 
precise ephemeris are insignificant. 
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homogeneous sphere, a "cannonball", and then trying to compensate for the 
obvious deficiencies of this model by adjusting, together with the Initial 
position and velocity, the coefficients that scale the drag and the solar 
pressure; this is usually done at intervals shorter than the total length 
of the arc, such as once per day (Marsh and Williamson, lb., 1980). A 
study of the consequences of using various simplistic, but practical, ways 
of handling the surface forces can be found in (Colquitt et al,, 1980). 

If A^ is the (pre-processed) altimeter measurement, h 8C the computed 

satellite ellipsoidal height, and Ah the error in the latter, then the 

s 

calculated value h wc of h w obtained from h c and A m is, according to (3), 

h - h - A 
wc sc m 

•= (h + Ah ) - (A + AA) 
s s 

= h + Ar - AA . (4) 

• w 

where the error Ar in the computed radial geocentric distance to the 
satellite has been put ip place of the ellipsoidal Ah , both being almost 
identical, to simplify the mathematical treatment. Replacing N and 
according to (2a) and (2b) in (1) and the resulting expression for h w in 
(4) leads to 

h - ((N - AN) + w + - AC- + w.) + Ar - AA . (5) 

WC C O be b t 


n 


Let 


Ah 


w 


"wo - (N c + W 


( 6 ) 


be the residual sea height * Subtracting N c and ? Gc from both sides in (5) 
one gets, after rearranging terms* 

Ah = Ar - AN - A 5 + w + w - AA , (7) 

W Got 

where A?_ = - (C + <5 + 6 '), 6 ' being the yielding of the solid Earth 

G Gc o 

to the ocean tide and 6 the solid Earth tide. 

On present evidence, and according to the way in which they have been 
defined above, AN, Ar, A? and w seem to have mostly long- wavelength 
spatial features (larger than 1000 km). The r.m.s. for each term in (7) 
is thought to be: 2 m or less (today) for both Ar and AN; nearly 1 m for 
w Q and for the geocentric tide; less than 0.1 m for w t (except over a 
small percentage of the sea surface); and 0.1 m (or less) for AA, provided 
there is a microwave radiometer onboard, as in SEASAT, that can measure 
the amount of water vapor along the altimeter beam, to estimate accurately 
then delay due to 'Vet" tropospheric refraction. 
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3.0 ORBITS 


3. 1 Orbital Motion and Keplerlan Elements 

As the orbits considered here are approximately elliptical, it is 
convenient to formulate their theory in coordinates that make the 
description of elliptical motion particularly simple. Among several 
possible choices, the most common is that of the six quantities known as 
the osculating Keplerian elements , which can be translated into the three 
Cartesian coordinates (x, y, z) for position and the three (x, y, z) for 
velocity, and vice versa, in an unambiguous way (see Table l(a-b) and Fig. 
1), Keplerian elements are used in the classical analytical theory of 
satellite geodesy. This theory has its origins in the old method for the 
study of planetary perturbations known as "Variation of Constants" (see 
Brouwer and Clemence, 1961), adapted by geodesists, in the early days of 
the "Space Age" (late 'fifties and early 'sixties), to the special task of 
mapping the complex gravitational field of the Earth. 

Before defining the "elliptical" coordinates, consider a Cartesian 
system with the origin at the center of mass of the Earth, or geocenter, 
the x and y axes on the plane occupied by the equator at the start of the 
orbit, or equator of epoch , the z axis aligned with the Earth's spin axis 
at that epoch, and the x axis lying along the intersection of the equator 
with the plane of the ecliptic, and pointing towards the vernal equinox. 

| The orientation is fixed with respect to the distant stars; in general, 


I 
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because the whole Earth, Including Its center of mass, is accelerated by 
external gravitational forces (the pull of the Sun, the Moon and the other 
planets), the system is quasi-inertial. The analytical theory is 
developed as if it were truly inertial, by making use of the idea of tidal 
potential as explained in the section "Tides ,, . 

The osculating ellipse is the two-body orbit that the satellite will 
begin to follow, from the point of view of an observer fixed to the system 
just defined, if, while driven only by their mutual gravitational 
attraction, both the Earth and the spacecraft were to shrink suddenly at 
time t, becoming point-like particles situated at their original centers 
of mass, but each retaining its velocity and momentum. The mass-center of 
the satellite would still have the same position vector _r with respect to 
the geocenter (the modulus |jr| being the geocentric distance r) and the 
same velocity vector f , which means that the osculating ellipse is tangent 
to the true orbit at time t (hence its name). Because of the nature of 
Newtonian physics (which is the one used here) given r_ and £ at t, the 
trajectory after t is entirely determined, so the osculating ellipse is 
defined by these two vectors. Furthermore, it lies entirely in the 
same plane as them, the instantaneous orbital plane , and has one 
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focus at the geocenter. The ascending node is where the satellite 


would cross the equator going north if It were to continue along the 
osculating ellipse, and the descending node is the point where it would 
cross that plane going south. The line of nodes , determined by these 
two points, is the trace of the instantaneous orbital plane on the 
equator. The closest point to the geocenter is the perigee, the 

furthest is the apogee ; they lie at opposite ends of the major axis of 

the ellipse, or line of apsides. The osculating Keplerian elements 
a, e, 1, ft, oj and M are: the semi-major axis a, the eccentricity 

e (0 < e < 1), the inclination I (angle between the equator and the 
orbit plane), the argument of the node fl (angle between the x axis and 
the line of nodes as shown in Figure 1), the argument of perigee (angle 
between the line of nodes and the major axis), and the mean anomaly M. 
The latter is related to the eccentric anomaly E (shown in Figure 1) 
through Kepler '8 equation 


M » E - e sin E , (8) 


^The reason why the osculating ellipse of the satellite has a focus at 
the geocenter is not the much smaller mass of the satellite, compared to 
that of the Earth. Even if both bodies had the same mass, or the system 
were attached to the center of mass of the satellite, rather than to the 
Earth's, the focus would still be at the origin of coordinates. This is 
because the osculating ellipse describes the instantaneous relative motion 
of one center of mass with respect to the other , whatever the masses 
involved, as long as the relative velocity is below the escape value (see, 
for example, Brouwer and Clemence (ib. , Ch. I, I960). An example is the 
"orbit" of the Sun in classical tidal theory (see paragraph "The tidal 
forces") • 
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a transcendental equation In E which must be solved Iteratively, given 
M. The mean anomaly would be the angular distance from perigee of a 
satellite moving at a constant angular velocity equal to n, the mean 
motion , 

, I / 2 

n * (pa ; (9) 



where ^ M(t) is the rate of change of M in the osculating ellipse , p is 

the product of G (the universal constant of gravitation) times the sum (M e 

+ M g ) of the masses of the Earth and the satellite. As M e » M g , this 

reduces to u = G M . If the Earth and the spacecraft were truly 

e 

particles, (9) would be Kepler's Third Law, and the orbital period would 
be T q = 2it n _1 . The ground track of a satellite is the line described on 
the Earth's surface by the point directly below the spacecraft, or 
subsatellite point. Table lc gives the relationship between the elements 
and <j> and X at the subsatellite point. Because of the rotation of the 
Earth and the motion of the satellite, the ground track is a (nearly) 
spherical helix wound up between the parallels of latitude * I. The pitch 
of this helix depends on the ratio between the terrestrial and the orbital 
angular frequencies. By tuning the orbit carefully, the helix can be made 
to close upon itself after a given number of days, forming a periodically 
repeating ground track which is also rotationally symmetrical with respect 
to the Earth's axis. This is done sometimes in Earth-surveying missions, 
including altimeter ones. 
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TABLE 1(a) 


CONVERSION FROM INERTIAL TO ELLIPTIC COORDINATES 


Given 


£ “ *£* + y£* + zR* 

£ " xi* + yj* 4* zk* 


where i*, J*, k* are unit vectors In the directions x, y, z. 


1. Compute a. Let 


and 


|r| « (x Z + y + z Z ) 
|*| - (x 2 + <f 2 + * 2 ) 


1/2 


1/2 


Then 


r P 


2 

2p-v r 


2. Compute e 


e cos E = 1 

a 


e sin E » _r • r / /pa 

/~ 2 2 ' ' ' 2 “ K 

e “ /e cos E + e sin L 


2. Compute M 


_ -1 ,e sin E. 

E * tg (— -) 

6 e cos E 

M * E - e sin E 


(E is in the quadrant where e sin E 
and e cos E have their signs as In 
the previous step.) 

(M > rr if E > ir) 


4. Compute 0 


-1 yz-zy 

a - tg 1 (V~r) 


xz-zx 


(Same comment as for E, regarding 
the signs of numerator and 
denominator inside the brackets.) 


5. Compute I 


_l i -zx+xz 
I - tg (• - -■■ — — - — 7 ) (I is always less than ir.) 


cos ft 


xy-yx 


6 > Compute u) 


« » tg [ 


(sinl)(x cosft + y sinft) 


] 


o *"1 / /l~*~ e _ E v 
- 2 tg (/-jl^ tan ^ 


(See comments for 
E and for 0.) 
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TABLE 1(b) 


fi 






CONVERSION FROM ELLIPTIC TO INERTIAL COORDINATES 

1. Compute the constants 

P x ■ cos 0 cos oj - sin to cos I sin ft 

Q x “ -cos ft sin to - cos to cos I sin ft 

Py ■ sin ft cos to + sin to cos I cos ft 

Qy * -sin ft sin to + cos to cos I cos ft 

P„ * sin I sin to 

z 

Q z = sin I cos to 

2 . Compute E 

Solve Kepler's equation 
E * M - e sin E 

iteratively, starting with E(0)« M. 


3. Compute xyzxyz 


r_ = £ a (cos E-e) + £ a / 1-e sin E 
r ■ — [q / 1-e^ cos E - p sin EJ 


where 


p = P i* + P j* + P k* and q = Q i* + Q i* + Q k* 

x “ y - z - *■ x "■ y z 


19 


TABLE 1(c) 


CONVERSION FROM ELLIPTIC TO SPHERICAL COORDINATES 

1 . Compute r 

r = a(l~e cos E) where 

3 

E - M + ( e - J-) sin M + ~ e 2 sin 2M + 0(e 3 ) 

(or solve Kepler's equation iteratively, as suggested in Table lb). 

2. Compute the true anomaly f (shown in Figure 1) 

f * M + (2e - j e 3 ) sin M + | e 2 sin 2M + 0(e 3 ) 

which is the equation of the center (for rough calculations, f^M is 
acceptable if e is small). 

3. Compute 

<j> = sin *[sin I sin(o) + f)] 

where sin - *(0) = 0 and <j> is in the same quadrant as OJff if 

7t IT 

- — < uH-f < j , otherwise it is in the same quadrant as ir-(uH-f). 

4. Compute A 

X - sin * [cos I sin (oH-f )/cos<t>] + 11-0 , 

where 6 is the sidereal angle of Greenwich, and A-JJ+0 is in 

the same quadrant as (ort-f) sign [■- - I] if I f — ; when 
I = , X-J2+0=O if urt-f < 7r , and X-JH-0 = ir if wff > n (of 

course, X is not defined at the poles). 

These expressions for <f> and X are the equations of the ground track . 


Knowledge of a, e, I, ft and w determines the size, shape and 
orientation of the osculating ellipse relative to the x, y, z axes* while 
M defines the position of the satellite in that ellipse and (together with 
a* e and p) also its velocity. Consequently, these six elements are 
equivalent to x,y,z, x,y,z, and can be used instead of them. The 
advantage of doing so is that they do not vary nearly as much along the 
orbit, with the exception of M, as the Cartesian coordinates and their 
derivatives do. There are other groups of six orbital variables used as 
coordinates (see Giacaglia, 1977, for example), but the Keplerian elements 
are sufficient for the purpose of this study. 

Satellites are sizeable, complex objects, so a full theory of their 
movements requires many more than six state variables. They spin, tumble 
(if not properly stabilized) and even change shape. Their detailed 
dynamics can be quite difficult to describe and doing so is the purpose of 
attitude control theory (see, for example, Kaplan (1976)). Any spacecraft 
considered here is designed to be sufficiently rigid so, once its attitude 
has been determined with the help of on-board sensors, the location of its 
center of mass relative to the altimeter is also well known. In such a 
"rigid" satellite, the position and velocity of the center of mass vary in 
time much like those of a material point with the same mass and driven by 
the same forces. For this reason, satellites are discussed here as if 
they were particles , and expressions such as "state", "position" and 
"velocity" actually refer to their centers of mass. 
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3.2 The Mean Ellipse 


Averaging the osculating elements up to time t over many revolutions 
of the orbiting body gives the mean elements a, <T, T, w and M at this 
time. These running averages describe the slow trends in the evolution of 
the orbit. They define a gradually changing mean ellipse . Usually, the 
main variations in this ellipse are: a steady, slow turning of its plane 
about the Earth's axis (precession of the line of nodes), a rotation of 
the major axis in the mean orbit plane about the geocenter (precession of 
the argument of perigee), and a slight departure from Kepler's Third Law 
(expression (9)) in the orbital frequency. ( 1) These changes occur because 
the Earth is not a homogeneous sphere, but a rather flattened ellipsoid. 

For an Earth-orbiting satellite, and except for equatorial orbits 
(I = 0) where the line of nodes is not defined, the precession of the 
nodes is given by the approximate law 

^ “ 2 " C 200 (if (1 “® } cosT » (10) 

_3 1/2 

where n * (p/a ) , R is the mean equatorial Earth radius, and C 2 QQ is 
the second zonal potential coefficient (this notation shall be explained 

(l^A further explanation of the behaviour of the mean ellipse is given 
later, in the paragraphs "The equations of motion in Keplerian elements", 
"The 'frozen' , repeat orbit", et seq. 


22 


later, see formula (26)).^^ Expression (10) shows that 0 is ruled mostly 
by the gravitational pull Qn the satellite of the equatorial bulge 
represented by CgQQ-* This pull adds a force normal to the orbital plane 
(except when the orbit is polar), causing this plane to precess much as a 
gyroscope does, and for the same reasons: to maintain the angular 

momentum vector constant in inertial space. Therefore, 0 f 0 unless 
I ■ j . For some values of a, e and I, ft ~ l°/day and the orbital plane 
completes a full revolution in one sidereal year. Such an orbit is known 
as heliosynehronous because it follows the Sun, keeping the angle between 
its plane and the direction Earth-Sun approximately fixed. A 
heliosynehronous orbit is chosen, sometimes, so that some onboard devices 
may maintain a proper orientation relative to the Sun. Possible examples 
could be the orbits of ERS-1 (Dow and Klinkrad, 1982) and POSEIDON (see 
POSEIDON report, 1983). 


The story for to is somewhat more complicated. If e is large enough, 
to is also governed by C 200 according to the expression 


2 —2 

to ~ j n C 20Q (1-e ) 11-5 cos 1 2 I] 


(ID 


so at I = arc cos yj\ ~ 63.4°, (or = 116.6°), known as the critical 
inclinations , to « 0 , changing sign as I goes through this value. 
However, for very low eccentricities (as those of the orbits of altimeter 


(1) R~ 6378.2 km; U ~ 398600.3 kra 3 /sec 2 ; C 2oo ~ -1082.63xl0" 6 ; the 3rd 

zonal C 3oo ~ 2.5xl0“ 6 , and C noo < C 3oo for all n > 3. 
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satellites) and depending on a, T and "e, the perigee may not precess at 
all but librate (oscillate slowly) about a mean value of tt/2 , or even stay 
for a very long time at ir/2 with no significant change, as in the 
so-called "frozen orbits" (see Cook, 1966, and Hough, 1981), which are 
important for altimetry because they can have very precisely repeating 
ground tracks (Cutting et al. , 1980; Dow and Klinkrad, ib. 1982). In the 
precessing case, the period 2 rr/ fib of a complete revolution of the perigee 
is known as the apsidal period . 

The oblateness of the Earth also affects M, making it slightly 
different from the Keplerian mean motion n: 

. o / vo 2 —3/ 2 „ 

M — n {l C 2Q0 (|j ( 1-e ) [3 cos Z X-lJ } . (9') 

Values of ft and to of ± 3° per day, and of M of 13-14 revolutions per day, 

are typical for the satellites discussed here, For a near circular orbit 
— -3 

(e^lO ), the orbital frequency (of successive passages through the 
ascending node) fluctuates slightly in each apsidal period about ft + to . 
Itself quite close to M (as to«M) , the angular frequency M+to appears in 
many important analytical expressions of orbit theory that will be seen 
later on. 

The combined influence of C 300 and other zonals makes orbits 

with e = 0 unstable , forcing them to become gradually more eccentric 

— -3 

(until e~10 for altitudes of about 1000 km, see expression (49)). 

2U 


Otherwise, changes in "a, 1?, and t are negligible, specially over the week 
or so that is the normal duration of the orbit arcs considered here. 
Variations in T, always very small, are somewhat more pronounced in the 
heliosynchronous case; they are caused mostly by the Sun's gravitational 
attraction. 

NOTE : While the time derivatives of M, "u, etc. are ft, w, etc., those 

of the osculating elements are M(t), *o(t), etc. 

An orbit is shaped not only by gravitation, but also by surface forces 
Such as air-drag and electromagnetic radiation pressure, which tend to 
change considerably the total mean energy E (kinetic + potential). With 
the sign convention adopted in physics, for spacecraft moving below the 
scape velocity (at which "a turns infinite) this energy is 

I = -p/2a . (12) 

Clearly, the decrease of E implies the decrease of 7T (notice that E is 
negative), which, in turn, means the increase of the mean velocity "v (for 
nearly circular orbits) because 
1/2 

v ~ (w/a) . (13) 

Therefore, dissipative forces like air drag tend to bring a spacecraft 
down while accelerating it. This is true as long as the braking force 
is gentle and the orbit remains close to circular, but is no longer the 
case when a satellite re-enters the denser layers of the atmosphere. 
The effect of drag on orbital motion has been explained in detail by 
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King-Hele (1964). As for electromagnetic radiation pressure (principally 
caused by the light of the Sun, either received by the spacecraft 
directly, or else reflected, or re-radiated as heat, from the Earth), its 
effect can be either a loss or a gain of energy, depending on the 
direction of the incoming rays, which may either push forward or brake the 
satellite. A paper by Rubincam (1982) gives a thorough description of 
non-gravitatlonal forces, besides drag and radiation pressure, and also 
lists many references on this subject. 

3,3 The Reference Orbit 

Altimeter satellites are placed in orbits of small eccentricity to 
keep the distance to the surface below always close to the optimum range 
of their instruments. They are also put high enough to make air drag 
small, but sufficiently low to have short orbital periods and, 
sometimes, finely spaced ground tracks. The Inclinations are chosen so 
most of the ocean surface is scanned. These and other often conflicting 
requirements result in compromise orbits that are, typically, about 

_ -3 

1000 km high (a ~ 7300 km), nearly circular (e ~ 10 ), and with 

inclinations larger than 60°. The orbital periods are all close to 
1.7 hours, or some 14 revolutions per day. Departure from circularity 
is of the order of ± 10 km. The ephemerides are calculated, usually, in 
arcs of up to one week in duration; their position and velocity errors 
nowadays are, according to available evidence, of the order of a few 
meters and a few millimeters per second, respectively. These errors are 
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sufficiently small to be studied using a first order orbit perturbation 
theory. The one introduced in this section is atuned to the particular 
characteristics of the orbits of altimeter satellites. It differs from 
the more conventional formulation (Kaula, 1966) in the choice of the 
reference orbit based on which the equations of motion and other 
relationships are linearized. 

In the standard treatment, the reference orbit is a precessing 

ellipse of constant e, a, and I with M, ft and w varying according to 

expressions (9')* (10) and (11). This orbit is governed by the Earth's 

’’central force term 1 * p/r and the attraction of the equatorial bulge 

alone. In what follows, instead, the reference orbit consists of a 

precessing ellipse defined by the arc averages a, e, I, fl, w and M of 

the corresponding computed osculating elements and their rates, together 

with starting values 0 , m and M at time t , the beginning of the arc, 

o o o o 

that give the best fits to 8, to and M, in the sense 
that oM- fi) (t-t Q ) ~ m(t), etc. The rate ft, in particular, is always quite 
close to the value given by expression (10). In general, unless the orbit 
is stabilized by on-board rocket engines, the mean elements and their 
rates will be different from one arc to the next, but this change should 
be very small. 

To see if the choice of the weekly reference orbit is reasonable, one 
can look at the actual orbit of SEASAT according to the ephemeris 
published by the Jet Propulsion Laboratory (Bauer, 1978). During its last 


month In operation (mid-September to mid-October of 1978), the satellite 

was kept in a ’’frozen", repeat orbit (i.e., with a repeating ground track 

and S) ~ 0). Over that period, a‘ decreased some 20 m per week, or 3 parts 

per million of its average value of 7169 km, while e" and T changed, in 

-3 

any of the four weeks, by less than 1 % of their monthly averages of 10 
and 108°, respectively. Likewise, !o fluctuated slightly each week about 
a mean of less than 1° per month; also the orbital frequency 
(approximately M + w ~ ft) , as indicated by the log of equatorial 
crossings, departed by less than. I part in six thousand from its overall 
monthly average of about one cycle every 101 minutes, while ft was close 
to 15° per week, and nearly constant. Therefore, the type of reference 
orbit adopted here can be a reasonable approximation to the true mean 
ellipse, over a weekly arc, at least for satellites like SEASAT. 

The radial orbital error is, at present, of the order of a few meters 
for altimeter satellites. The mean ellipse, on the other hand, can be 
several kilometers away from the true position of the spacecraft, mostly 
because of short period perturbations due to the large second and third 
zonals (caused by the oblateness and "pear-shape" of the Earth). Is it 

valid to describe, at least qualitatively, the radial error using a first 
order perturbation theory based on such a reference orbit? 

To separate the useful signal from the radial orbit error At it is 
more important to know well the frequencies and phases of the spectral 
lines of this error than the precise values of the amplitudes . From this 
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point of view, the choice of reference orbit made here is probably 
adequate. In fact (as shown later), the corresponding perturbation theory 
seems to explain quite well the results of some computer simulations, at 
least qualitatively. 

3.4 The Radial Orbit Error 

The error that matters here is that in the geocentric distance r, 
l.e. , Ar in expression (7). In osculating Keplerian elements, this 
distance is 

r » a(l - e cos E) , 

where, expanding E as a function of M, as in Smart (Ch. V, par. 71, 1931 
(6th ed., 1977)), one gets, to the order of (or 0(e^)), 

cos E = (1 - 4 e^) cos M - ^ e cos 2M + e^ cos 3M 

so 

v - a(l - e cos M) + O(e^) . (14) 
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Because of the small eccentricity of an altimeter satellite (~10 ), the 

second term In (14) can be neglected, and 


MwtJLitadfc- -Mh 

i 


i 

1 


r » a(l - e cos M) . ( 15) 

This equation applies also to the mean ellipse, so the maximum radial 
departure of the mean orbit from a circle is ae, or some 10 km 
for e~10 _ "* and a~lo\m. Differentiating (15) gives the first order 

approximation to the radial error 

(16) 


(17) 

and t" - is the time (or epoch) when the k weekly arc starts; the overbars 
o 

indicate that the partial derivatives are taken on the reference orbit 
jj ■ ■ corresponding to this arc. 

ij ■ ; : , •: : ■ / ' : ■ . ' ■ : ; ^ ■ ’ v / ; : : : w : V 

ft 


Ar = Aa - (Aa e + Ae a) cos M + AM a e sin "M , 


where 


TJ = M(t - t£) + M(t£) 

j . ••• 

i- ■' ; 'V: : ; " ; ; . ■ v v ' 
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3.5 Temporal Characteristics of the Error 

The frequency spectra of A a, Ae and AM can be separated Into four 
parts. The first three are sets of well defined lines with frequencies 
ranging: (a) from zero to a small fraction of a cycle per year; (b) from 

a small fraction of a cycle per year to a cycle every few months; (c) from 
a cycle every few months to about 50 cycles per day. Superimposed on 
these lines there is (d) a continuous background, due mostly to the 
inaccuracies in the modelling of the surface forces. 

(a) Practically Constant Part : Over one week or less, any 

oscillation in the lowest range of frequencies is indistinguishable from a 

constant, and it will be regarded here as such, if Aa, Ae and AM in (16) 

are constants, the resulting Ar consists of an offset between the computed 

and the true orbit, equal to Aa, and of an oscill&tion at frequency 

A, which is close to (M+u), or about one cycle per revolution. The 

_ _ o . 9 

amplitude of this oscillation is [(Aae + Aea) + (AMae) ] . 

(b) Nearly Secular Part : Oscillations with periods in the second 

band may show appreciable variations during one week, but their rate of 
change will vary so slowly that they will resemble secular changes. If 

• j£ • • 

t Q is the time at the beginning of the k weekly arc, these changes can 

jj. : . 2 

be approximated very well by expressions of the form c (t-t )+d (t-t ) 

o o 

(where cf* and d^ are constant and c^ « d^) for all the Keplerian 
elements. Both c^ and d^ will vary slightly from one arc to the next. 
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The plots of the corresponding Aa, Ae, AM, etc#, versus time will depart 
from straight linos just enough to show that they are not truly secular. 
These slow variations must be reflected in Ar, according to (16), as a 
slow change plus an increasing oscillation of frequency ft with a slightly 
bent (quadratic) envelope. As w is usually very small, ft will be very 
close to one cycle per revolution. 

(c) Periodical Part : Mostly caused by errors in the geopote ( r.tial 
field model, it consists of sinewaves whose frequencies depend on the 
rates h, w and ft of the reference orbit, as well as on the spin rate of 
the Earth, 6, Because these frequencies are higher than those of the 
errors mentioned previously, the oscillations in this band will be called, 
in what follows, "shorter period" errors. Other periodic errors are 
produced by uncertainties in the tidal models; these are smaller and have 
different frequencies, in general, than those due to the field model. 

(d) Other Temporal Variations : While they may also include secular 
variations and oscillations, errors caused by inadequately modelled 
surface forces, etc. , do not have sharply defined lines in their frequency 
spectra. During one week, however, probably much of their total effect 
on Ar not included in (a), (b), or (c) can be represented by a low degree 
polynomial (a cubic or a quartic) of the general form: 
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( 18 ) 





(ng) 

Jk 




where J Is a small Integer, (ng) denotes "non~gravitational, M and 
the r^®^ are constants which may vary from arc to arc ( "k" is' the 
sequential number of the arc). 

Errors in the calculated precession and nutation of the Earth modify 
some of the lines in the spectrum of hr related to the diurnal tides. The 
periods of any significant effects are long enough for them to be lumped 
together with the constant, secular and polynomial parts (a), (b) and (d), 
and, in any case, their amplitudes are very small. Therefore, they will | 

not be considered here further. For details, see Reigber (1981) (he j 

considers the total effects of precession-nutation on the actual orbit), } 

Having outlined the features of the radial error, it is time to ] 

consider the mechanisms that shape them. 



3.6 The Equations of Motion and Their Linearization 


Let V be the part of the gravitational potential off the Earth (or 
geopotential) which is Invariant in an Earth-fixed system of coordinates 
(the time-varying part due to the tides will be discussed in "Tides"); the 
gradient VV is the corresponding gravitational, acceleration. If is a 
vector representing both tidal and non-gravitational accelerations, then 
the Newtonian equations for the motion of the center of mass of a 
satellite in an inertial frame can be written in vector form as 

• • 

tm VV (r, t) + f ; (19) 

V is given here as a function of both position and time, because the field 
is rotating along with the Earth. One interesting consequence of this 
time dependence is that the total energy (kinetic+potenfcial) of the 
satellite in an inertial system of coordinates is not conserved, in 
general . Only a central force field, or a zonal field, both off them 
invariant with respect to a rotation about the z axis, are conservative; 
not so the more irregular field of our planet, except in Earth-fixed 


coordinates. For a discussion of this problem, see Hotine (Ch. , 28, 
1969). In general, f_ is also a function of _r , t, and r. 


The accelerations in (19) are nonlinear functions of _r_ (think, for 
example, of the simple case of the field of a particle, 
where |w| ® pr ). For the treatment that follows, (19) has to be 
expressed in terms of the Keplerian elements and of 0, the sidereal angle 
of Greenwich that accounts for the Earth's rotation. Lagrange obtained 

his transformation of the equations of motion by a very ingenious and 
rather laborious process described in Brouwer and Clemence (ib. , 1961) and 
also in Kaula (ib., 1966). The end-product is a system of six first order 
differential equations, Lagrange's Planetary Equatiohs, one for each 
osculating element "s^” (i = 1, 2... 6); as used in geodesy, they are of 
the general form 


s « L (a, e, I, ft, a), 9, M, 6, s k ) , (20) 

IS, " o 


^^The potential part of the total energy depends, like V, on the position 
of the satellite relative to each of the particles that form the Earth. 
Such relative positions and the resulting V are coordinate invariant . 
The kinetic energy , on the other hand, is proportional to the square of 
the velocity, which depends on how the spacecraft moves with respect to 
the reference frame, so it is not coordinate invariant . In any reference 
frame, however, the sum of the potential and kinetic energies of the 
spacecraft and the Earth are constant. To keep the books straight, the 
Earth must change its own velocity in the chosen frame, to compensate for 
the variations in the kinetic energy of the satellite; our planet being, 
by far, the more massive of the two bodies, such changes are 

imperceptible. The only system where an orbit must have always a constant 
total energy is, therefore, an Earth-fixed system which moves and rotates 
with the planet, so the velocity of the planet is always zero. 


where 

L 

8 i 

is a nonlinear function, 3 is 

a 

vector 

of 

parameters 

that 

1 

appear 

in the expressions of V and fj 

and 

k . 
s is 

~o 

the 

vector of 

the 

initial 

k 

conditions at t * t (the 

Btart of 

the 

k arc) : 
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* (a^ e^, 1^, , M* 1 , ^^-6^). These equations are shown in detail 

-O v O* O* o’ 0*0*00 1 

later (expressions (36)). 


To first order, the error As^ in s^ 

k k 

As . in the components 3 . of 3 and s , 
io r j — io 

approximated by the differential of which 


due to small errors 
K 

of s q , respectively, 

, according to (20), is 



can 


and 

be 


As 


s i - \ V <V 46 j + \ D s (L S > 4s „o • 

j j i u=l uo i 


(21a) 


where the derivatives are taken at and 


3L . 3L 3s 

s. 6 s u 

V <V = *33T~ + \ ~Ts~ W~ 

j i j u=l u j 


2ii 

* 


(21b) 
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while D uo (L s ^) has a similar expression. Equation (21a) Is the general 
form of the linearized equations of motion , and (21b), that of the 
variational (or sensitivity) equations. The variationals, which make up 
the "skeleton" of the linearized equations, are solved often by themselves 
to find the values of the unknown 3s / 3fi , . These values, in turn, are 
used to calculate the elements of the matrix of "partials" needed for 
orbital adjustments, force model improvements, sensitivity studies, etc. 
The sura in (21b) represents the dynamic terms, and the forcing 

terms of the variationals. 


3.7 The Homogeneous Linearized Equations 

The differential equations given by (21) are thoroughly coupled to 

9L s 

each other through their "interaction terms" -gg— ^ , which multiply the 
3s u u 

unknown As g = -gzr- A3-j . They are also time-dependent, in general, like 
3L p j p j J 

the -gg — themselves. These characteristics make the exact solution 
possible only by numerical integration with electronic computers. To 
obtain the analytical expressions needed for the present study it is 
necessary to resort to some simplifications, This limits the numerical 

accuracy of the solutions, but leaves virtually intact their qualitative 
properties, which are the ones of real interest here. The usual approach 
is to ignore the disturbing forces represented by £_ in (19) (including all 
dissipative forces like drag) and also the departure of the terrestrial 
field from that of a perfect sphere (or of a point-like mass), thus 
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disregarding Earth-rotation as well. This dissipation-free, spherical 
approximation , can be accurate enough to describe the main characteristics 
of the errors in precise ephemeris calculated with force models that take 
already into account most of f and of the "non-spherical" part of V, so 
the main long period effects (particularly regarding A and &) are 
included in trie mean ellipse of the computed arc, which is also the 
reference orbit here. The simplifications are more admissible when the 
inclination I is high and the eccentricity e is small, as in the case of 
altimeter satellite orbits, because most of the main terms ignored are 
proportional to e or to cos I (see Kaula, lb. , 1966, Ch. 4, par. 4.3). 


In a spherical, or central, force field, when _f = _0, the orbit of a 

satellite would be a simple ellipse obeying Kepler's Laws ( two-body 

orbit). The elements a, e, I, m and would be constant for all time, and 

3 1/2 

only M would change according to the Third Law: M = n = (y/a ) , where 

"n" is the mean motion. Consider a "true" elliptical orbit and a 

"computed version" of it obtained With the exact value of u, (i.e. , the 

field is known), and assume that they differ slightly in the initial 

values of their elements. Then AM(t) = M( coraputed ) - M (true ) 

3M 3 n 

= — Aa(t) = - — — Aa(t). Since none of the other s. depend on 

9a 2 a i 

t, As =0 for all of them. As now the osculating ellipse must be the 

. „ - 1/2 

same as the mean ellipse, a = a and n = n = (y/a J ) . Therefore, the 
effects of the initial state errors on the computed a(t) , e(t) and M(t) 
must satisfy . 
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Aa(t) * 0 


(21 'a) 


Ae(t) - 0 (21 'b) 

AM(t) - - 1 JE1 Aa(t) , (21 'c) 

2 a 


These are three of the six homogeneous , or unforced , linearized equations 

of motion (i.e., all 3L_ /38. *» 0 in (21b)) for the spherical, 

s i J 

dissipationless case. For I, w and J2 the corresponding unforced equations 
are all of the form As^= 0, like (21'a-b), but they are not relevant to 
the study of Ar (see expression (16)) , As the linearization is made along 
the computed orbit, both "a and "n correspond to this orbit. 


3.8 Orbit Error Due to the Estimated Initial State 


Each arc of the computed orbit is fitted to the tracking data by a 
least squares adjustment of the initial conditions(^ (the s^, where "k" 
is the number of the arc) and of a few other parameters, such as drag and 


^^"As all our observations, on account of the imperfection of the 
instruments and of the senses, are only approximations to the truth, an 
orbit based only on the six absolutely necessary data may still be liable 
to considerable errors. In order to diminish these as much as possible, 
and thus reach the greatest precision attainable, no other method will be 
given except to accumulate the greatest number of the most perfect 
observations, and to adjust the elements, not so as to satisfy this or 
that set of observations with absolute exactness, but so as to agree with 
all in the best possible manner." 

Carl Friedich Gauss, "Theoria Motus" (1857). 
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radiation pressure coefficients. The data, the coordinates of the 

tracking stations and the various force models used to integrate the orbit 

are not perfect, resulting in adjusted values that are wrong to some 

extent. To understand the radial errors that an incorrect initial state 

produces in the computed orbit, one can solve the approximate linearized 

homogeneous equations of motion (21') for Aa, Ae, and AM with the 
■ Ic 

errors As^ q as initial conditions? 

Aa(t) * Aa k (22a) 

o 

Ae(t) = Ae k (22b) 

o 

AM(t) = - 4 - A a k (t-t k ) + AM k , (22c) 

2 — o o o 

a 

k 

where t is the starting time of the arc. Kaula gives more complete 
o 

expressions (ib. , Ch.4, equations 4-25) by including small "interaction 
9L si 

terms" — approximately proportional to C^qo tllat have been neglected 
uo 

in (21'a-b) because they would be zero for a spherically symmetrical 

Earth. As shown in (22), the adjusted initial conditions, by being 

slig'htly wrong, can cancel some of the constant errors in a and e, and 

some of the constant and secular errors in M. Replacing Aa, Ae and AM 

(ic) 

from (22) in (16), the contribution Ar of the initial conditions to 

the total radial error is 


Ar(^ c )(t)*AaQ+[-(& ao * c e+A e ^a) C os(MQ-M tg) - ^ 5" Aa^ t^+AM^ae 

sin(fl£-Mt£)]cos Mt+[(Aa^eH-AeJa)sin(Hj-M t£)-(-| •§ Aa£ t^+AM^)ai 

cl 

cos(RQ-Mt^)]sin Mt-[ae Aajj t^sinCHg-Mt^Jt cos Mt 

a 

- [ae J 2- Aa£ t£cos(Rjj-Mt£)]t sin Mt 

. & ■ 

- Aao+A^ ic ^ cosMt+B^ lc ^ sinMt+C^ lc) tcosMt+pj^ tslnMt , 

(23) 

, ,(ic) _(ic) _(ic) , (ic) . . , c 

where A , ,B , ,C . and D . are the values of the square 

K iC K K 

brackets. In general, M will change slightly from arc to arc, due to 
orbital decay. 

3.9 The Free Response of the Linearized Equations and Resonance 

Expression (23) gives the radial component of the u nforced , 
homogeneous, or natural response of the linearized equations of motion. 
As the linearization has been done ignoring the term f_ that, in (19), 
accounts for all dissipative forces, (23) describes the behaviour of an 
undamped dynamic system, whose unforced response includes a non-decaying 
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oscillation of angular frequency M, the natural frequency . In reality, 
there is damping due to drag, etc., but for orbits higher than 200 km this 
damping is very light. For altimeter satellites, where e” ~ 0, the system 
is virtually a stable harmonic oscillator . ^ ^ The along track and across 
track free responses consist of similar oscillations, plus a slow secular 
drift along track, so the complete error due to initial state 
uncertainties makes the computed position ’'circle" the true one in an 
elliptical path once in every revolution of the spacecraft, while also 
moving slowly away from it along the orbit. When driven by a disturbing 
force that has a periodical component at a natural frequency, the forced 
response of an undamped system must contain growing oscillations of the 
same frequency. This is known as resonance . In the present case, this 
happens when the ephemeris crosses the same disturbance at repeated 
intervals equal to the period of M, or about once per revolution. As 
shown later, certain errors in W and in £ can produce this effect, 
which may grow into a large perturbation over a sufficiently long time. 

The resonant character of the linearized equations is not plain from 
their formulation in Keplerian elements, but is quite clear when they are 
given in terms of the perturbations of the radial , along track and across 
track components of the position vector, and of their first time- 

The "unstable" last two terms of (23) appear only in a first order 
approximation to 4r. The actual effect is an oscillation of frequency ft 
modulated by a periodic envelope whose fundamental frequency is 4ft, the 
error in M. As 4ft « M, over a sufficiently short interval of time (in 
which the perturbations remain small) the envelope seems to be expanding 
very slowly, at a steady rate. 
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derivatives. Written in this way, the equations are sometimes known as 
Hill's equations (see Kaplan (1976), and also Colombo (lb., 1984, Ch. 1)). 

The best known and most striking effect of resonance in the whole 
solar system is the series of gaps that divide the rings of Saturn. 
Inside each gap, any orbiting particles would have a period congruent with 
that of a major moon, but the gravitation of this moon disturbs greatly 

the motion of such particles, sweeping them out and keeping the gap open. 

/ 

/ 

* 

3.10 The Spherical Harmonic Expansion of the Geopotential in Keplerian 
Elements 

Outside any ideal geocentric sphere S that contains virtually all the 
matter of the Earth, including that in the atmosphere, V (the time- 
invariant part of the potential in Earth-fixed coordinates) can be treated 
as a harmonic function satisfying Laplace's equation 

V 2 V = 0 , (24) 

which expresses the conservation of the flux of lines of force in empty 
space. 
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Let Che spherical coordinates (r, <j>> A) correspond to the geocentric 
distance, latitude and longitude in an Earth-fixed equatorial 
system, A^O being the longitude of the Greenwich Meridian. A spherical 
harmonic of degree n and order m has the form 

Y nma ^* x ) * P nm ( sln cos (m A - — a) , (25) 


where P nm is the associated Legendre function of the first kind with the 

same degree and order. The Y (<j>,A) are orthogonal in the sense that the 

nttjot 

integral of Y (f r X) Y- „ * ^(<j>,X) over the unit sphere is zero if 

nmot n m ot . , 

R n+i 

n £ n', or m A m' or a f a'. The product (-) Y ( d> , A) , known as a 

solid spherical harmonic , is a solution of Laplace's equation, and any 
other solution, V for example, can be expanded in a series of solid 
spherical harmonics 



00 


1 

n~0 



m^O 



a=0 


c (-) 

nma x r 7 


n+1 


Y „(M) 

nma ST> 


(26) 


outside any Earth-enclosing sphere S (see Hobson, Ch. Ill and IV, 1931, 

and Heiskanen and Moritz, Ch. I and II, 1967). R is usually the mean 

equatorial radius of the Earth; C is the (dimensionless) spherical 

nmot 

harmonic potential coefficient of degree n and order m. By definition, 
Cqoo = 1 an d all C n01 =0, while all the C = 0 in a geocentric system 
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(r, 4>, A).^ The geoidal height N has a similar expansion, with both r 
and p/R replaced by R (this spherical approximation is quite adequate 
for the relatively low degree terms that dominate the geoid). The 
expansion of the geoid inside the sphere S (i.e. , on the reference 
ellipsoid) is valid, but the reason for this is subtle (see Moritz, 
1980). The error AV in V is 


AV 


M. y 

R *• 


n 

l 


l AC 


n=*2 m^O a®0 


nma 


<“) 


n*f 1 




(27) 


where the AC are the errors in the coefficients of the field model, and 
nma 

they are zero for 0 < n < 1. 

To study the motion of satellites, these formulae have to be 
converted to osculating Keplerian elements; the details of the 

transformation are given in (Kaula , ib. , 1966, Ch. 3). The variable 


Q' = n-e , 


(28) 


^ If the Cartesian - ? axis Were always aligned with the Earth's main axis 
of inertia, the C,^ = 0. In the actual instantaneous field they are not 
zero but very small and time-varying, depending strongly on the definition 
of the "Earth-fixed" equatorial system (see Reigber, ib. , 1981). Their 
mean values (which enter in the expansion of the time-invariant V) are 
also.- very small. 


where 6 is the Greenwich sidereal angle, must be used Instead of SI to 
take into account the rotation of the Earth. Introducing momentarily the 
true anomaly f (Figure 1), expression (25) can be written as^ 


Y „,n„<*>» * \ W l) l<n-2p> <“*» ' * 9 '> 

p»0 


l F <» 


p»0 


nmp nraapq=0 


(m,M, 0") 


(29) 


for small a, where 


c nma pq( u »M, 9'') “ cos {(n-2p+q)(Urt-M)~qaH-m0'*) - 5 ( af|(l-(-l) n m )J). 


(30) 


Expression (29) will be used later on to formulate the ocean tides in 
Keplerian coordinates. Solid spherical harmonics can be written as 
functions of the osculating elements of a satellite at a point (r,f,X) as 
follows (Kaula, 1961): 


D n+1 „ n+1 n °° 

& Y „„a« , - X) ■ <4> i l «.»>> c C»,M,9') , 


p=0 


nmp 


qja-e 


npq nmapq 


(31) 


d^The approximation above would become exact if e=0. From the "equation 
of the center" and the expressions for $ and X at the ground track (all 
given in Table lc, see also Smart, ib., 1931, ch. V), the maximum along 
track departure between a slightly elliptical orbit and a circular one of 
the same period is rfu ~ 2ae (consider the simple case where 1=0, in 
Table lc). For e ~ 10“^ and a — lO^km, <fu — 20 km, so (29) must hold well 
up to degree n — 200, l.e., up to spatial wavelengths of some 200 km. 
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The F nm p(I) and G n pq(e) are commonly known as inclination and eccentricity 

functions, for obvious reasons; the G (e) are also called "Hansen's 

coefficients". There are expressions for computing the and the 

nmp 

C 6 ) 

G , some of which can be found in (Kaula ib., 1966, or in Giacaglia, 
npq 

1977); they are not needed here. 

For near circular orbits (small e) , two important properties of the 

G(e) are: that they are approximately proportional to e”l c ll, so they 

npq 

decrease very fast with increasing |q|, and that 


him G n pq(e) _ |1 if q = 0 
e + 0 0 otherwise 


The expansion of V or AV, in Keplerian elements, takes the form 


V 

{or} 

AV 


“ n 1 „ n+1 

£ l l l <f> 

a r\ A « 


nma 


rO i m*0 cs»0 

n* to I 


{ or 
AC 


} l t 


(I) 


nma 


l G (e)c (w,M, 6'), 
p =0 nmp q*-» n P9 nmapq 


(32) 


Because of the altitude attenuation factor (R/a)^ n+ ^, the higher the 
degree of a spherical harmonic coefficient , or of its error, the lesser 
its effect on the potential and on the gravitational acceleration acting 
on the satellite. As a consequence, for computing near circular orbits, 
the series for V can be truncated roughly at degree n - 2xR/(a-R) , 


or n ~ 40 when a - R - 10 km. 


Also, except for a few coefficients, 
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mainly of low order (n < 10) or of "shallow resonant" order (n=13-*15, 

26-30. 39-45), the orbit errors due to the AC tend to vary in 

nmo ■ 

proportion to the attenuation factor (R/a) n+1 and the size of the AC 

nmo 

The latter ones, in today's models, appear to increase with n and m until 
they reach nearly 100% of the at n=m~30. Above this degree and 

order, the combined effect of all the estimated coefficients reproduce the 
data (tracking, altimetric ocean heights, gravimetry, etc.) reasonably 
well, but their individual values cannot be trusted. The size of the 
actual coefficients follows the approximate law 




^ (n+m) f 


1/2 

\ 


where 6 - is the delta Kronecker and 
m0 


„ i/2 


*£ < 2 n + l )‘/ 2 


(33) 


according to "Kaula's rule of thumb" (Kaula, ib. , 1966, Ch. 5, and Kaula, 

1967); o is the degree variance 
n — 



n 1 (n+m)! (C )^ 

f 2 nma 

m=0 ot=0 (n-m)!(2n+l)(2-6 A ) 

mu 


(34) 


and equals the mean square value (or "power") of the sum of all , harmonics 
of V of degree n over the whole of the Earth's surface (see Kaula, 
1967b). The quantity {a^/ ( 2n+-l ) } is approximately the average size of 
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the fully normalized coefficient C . (whose square is the general term of 
r - nma 

the sum in (34)). These coefficients correspond to the fully normalized 
harmonics, scaled versions of the Y , the integrals of whose squares on 
the unit sphere are all unity. The define the spherical harmonic power 
spectrum of V. "Kaula's rule", based on studies of terrestrial gravity 
measurements and of early field models obtained from satellite tracking 
data, has been shown to be a reasonably close guess of the actual power 
spectrum of the geopotential for n as high as 200, although its values for 
the o n are, on the whole, rather high between degrees 8 and 60, and rather 
low above 60. This conclusion is the result of a number of global 
analyses of altimetry and gravimetry. Wagner and Colombo (1979), and Rapp 
(1979) have given more accurate formulations for the o^, but "Kaula's 
rule" still has the convenience of its greater simplicity. 

3.11 The Equation of Motion in Keplerian Elements 

In Keplerian elements, the equations of motion of the satellite are 
known as Lagrange's planetary equations. If E is the total energy 

(kinetic + potential) at any given time, and F=-E, then it can be shown 
that 

F = JL- + v - ii , (35) 

la r 
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where F, r and a are Instantaneous values* The Lagranglan equations are 


. 2 3F 

" na 3M 


(36a) 


e(t) 


2 2 1/2 
(1-e ) 3F (1-e' ) 3F 


2 2 

na e 3M na e 


3 CD 


(36b) 


M(t) = - 


l-e^ 3F 2 3F 


na e 3e na 3a 


(36c) 


5>(t) 


cos I 3F (l-e) 1/2 _3F 

n « /n A * *» A n 


2,. 2.1/2 . T 31 T 2 3e 

na (1-e ) sin I na e 


(36d) 


i(t) 


cos I 


3F 


3F 


na 2 (l-e 2 ) 1/2 sin I 3o) na 2 ( 1-e 2 ) 1/2 sin I 3fl 


(36e) 


ftco “ 172 -rr > 

na (1-e) ' sin I 31 


(36f ) 


-3 l f 2 

where n = (p/a ) . 


The largest effects produced by the gravitational field come from the 
central force (Cqqq) and the oblateness (C200^* The first gives orbits 
their general elliptical shape, and the second causes most of the secular 
variations in their mean ellipses. Careful inspection of (32) shows that 
^200 assoc i ate< * with a term in the expansion in Keplerian elements of V 
that does not contain the fast changing M. This term 
is V' » (1 c|) 2 F 201 ( I > G 2io( e >- According to Kaula (ib. , 1966) : 
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*201<*> - 1 sl ° 2 1 - i 

and 


GzioCe) 


O-e 2 ) 


-3/2 



Replacing V with ~ Cqqq + V* in (35) and the resulting 
equations of motion become, taking the formulae for F 2 qj 
into account. 


a<o - o 

e(t) = 0 

i(t) = o 

In r riv 

2 C 2.00 V 

3 n /In 

4 n C 200 ( a ) 

^ fc ^ * R ~ "4 n C 200 


n(t) = 
5>(t) = 


COS I 

2.2 

(1-e ) 


[1 - 5 cos I] 


/Rv 2 [3 cos 2 I - 
V 3/2 

(1-e 2 ) 



F in (36), the 
(I) and G 2 ^Q(e) 
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The first three expressions show that the oblateness cannot cause secular 
changes in a, e or I. Replacing instantaneous values with mean values 
(which are very close to the former, in general) the last three 
expressions become (10), (11) and (9')i respectively. 

3.12 The Forced Linearized Equations 


Clearly, Lagrange's equations are not linear in their unknowns, the 

Keplerian elements. Linearizing according to (21), but assuming that the 

only errors are those in the C of the field model, sc 

nma * 

k 

all Asj = 0, produces equations of the form: 


, \ « n 1 

41 8> = l l l 

* rv 


3L 


3L 


n=2 m=0 a=0 


3C 


nma 


AC + 
nma 


6 

l 

j=l “s 


~ 4s j 8> 


j 


(37) 


where the superscript (g) indicates the gravitational origin of the 
perturbation or error (i.e. , the AC nma ). 


Notice that the sum with respect to "j" in (37) is the dynamic part 

of the linearized equations (21), while the (3L 0 /3C„ m/v ) AC „ are 

■ ' Sj nma' nma 

forcing terms. The linearized equations are mutually coupled by 
"interaction terms" of the form 3 S As j®- where the derivative can be a 
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function of time. This time-dependence makes an analytical or 
closed expression for the solution impossible in general. Fortunately, 
as pointed out when deriving the approximate homogeneous form of 
the equations, the terrestrial field is close to that of a spherically 
symmetrical body (a point mass, for example) where F * p/2a 
is constant and all interaction terms in (37) are zero except 
Aa^ g ^= -(3 n /2a) Aa^ g ^ in the equation for AM. Those vanishing 
terms are, in the actual field, of the order of C^qq* or about 10“^ times 
smaller than, than rest of (37), and it is quite reasonable to ignore them 
for the purposes off this study. 


To arrive at the equations for Aa, Ae, and AM, which are the 

ones needed to find that of Ar, replace V - p/r in (35) with its 
expansion in Keplerian elements according to (30) and (31) and then 
carry out the partial differentiations indicated by (36) and (37) 
at the reference orbit, where M=M + M(t-t Q ), + fi>(t-t 0 ) and 


3L C 


0" = n o -0 o + (n-0)(t-t o ), ignoring all 

The result is: 


As 


(g) 


3L, 


M aJs) 


except Aa 
9a 


2 

M(t) (g) = 2p/(n a ) 1 AC /^) n F (I) G . (e) (n-2p+q) 

nraapq nraa W ™P n Pq 


cos I ( (n-2p+q ) (a*M) - q ayfmS ')t + ^ cq)q£) + j] 


(38a) 
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/ v _ _3 , v \ n _ .. _2 1/2 2 1/2 

Ae(tr g; «u/(n a ) £ AC nma (-] F (I)G (e) [( 1-e ) (n-2p+q)-(n-2p)](l-e ) 

nmpqa ' ' v 


cos[((n-2p+q)(5H-ft)-qfiH-m6 , ')t+ij) + -^1 

nmapqo 2 


(38b) 


M(t) (g) « u/(n a^ n ^ a A c nma (f)^-( i | i -)(3Gnpq(i)/3e)+ 2 ( n +l)Gnpq<e)]Fnmp( I ) 


cos [ ((n-2p+q )( o>fM)-qa>fm0'')t+(|i 


l-2| i a(t) <S) , 


nmapqo 2 a 


where 


(38c) 


<}>„ m „* ( n " 2 p+q ) ( u +M ) - q to +m Q' 

nmapqo o o o o 


\ ir[a + { (l-(-l> n_m ) J 


9 ' = $ 1-9 


(40a) 


n = (p/a 3 ) 1 / 2 


(40b) 


M = M(t k ) - o)t k 
o o o 


(40c) 


and similarly for oj^ and 8*; t* is the starting time of the k arc. A 
different choice of starting time can be made in the case of "frozen", 
repeat orbits (as explained in the comments following equation (52)), to 


show, more clearly the main properties of both Ar and the tidal signal in 
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altimeter measurements taken along such orbits, which have precisely 
repeating ground tracks. 


Expressions (38a-c) are the forced or unhomogeneous linearized 

equations. Notice that if all the AC are zero, the (38) become 

nracx 

identical to the unforced or homogeneous equations (21'a-c). While the 

former are based on a precessing reference ellipse, as shown by 

M, a) and %* in the arguments of the cosines of their forcing terms, the 

(21') are not. However, both sets of equations agree perfectly with each 

other because both involve the same "spherical approximation" that 

neglects among their dynamic terms those proportional to C 2 qq linking 
(g) (g) (g) (g) (g) 

Aa , Ae , AM to Am , Aft and, thus, to any secular variations 
in to and ft whatsoever. The terra - j (n/a) Aa(t)^ is often ignored in 
the formulation of these equations, because the factor (n/a) is very 
small, so only very large perturbations in the semimajor axis would have 
an effect on AM^ , anil leaving it out makes the integration of the 
equations immediate. Here this term has been retained mostly for the sake 
of consistency. 


To solve (38a-c) the first two equations are solved by direct 


integration respect to time in the interval t < t < 


Aa (g) (t) = f. Aa (g) (t')dt' 

t* 

o 


Ae (g) (t) - Ae (g) (t')dt ' . 

t K - 

o 
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Next, Aa^(t) is replaced in (38c) and this is, then, integrated to 
obtain AM (g) <t). 

( & y 

Clearly, all solutions As^ e/ (t) obtained in this way satisfy the 
condition 




(41) 


This is physically meaningful, because, due to inertia, the finite change 
in the driving force due to the AC cannot affect the ephemeris 
instantly. 


3.13 Orbit Error Due to Inat':?\rect Potential Coefficients C 

nma 

When integrating the various terms in (38), there are two situations 
that must be considered separately. 


(a) (n-2p+q) (tirt-M) -qort-mO^ 0 


This is the normal case 


the corresponding part As. of the 


solution, for Aa and Ae is of the form 
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I 

nmapq(nonres) 


AC s cos[((n-2p+q)Cwfft)-qaH-mS > )t 
nraa inmpq 


+ <*> ] + As^ 

T nmapqo io 


(42) 

with 

s * s [(n-2p+q)(a>^M)-qu>^me•' , ] , (43) 

inmpq inmpq 

cvjr 

where the s^ nm pq are the functions of a, e and I in (38) and the As^ q 

satisfy the inertia condition (41). The symbol " \ " denotes 

nmapq(nonres) 

the sum of all terras where the frequency is not zero. As for AM^ S \ after 
replacing Aa^ in (38c) according to (42) and (43) and integrating, one 
gets 

4M<S> . iff- . 1 I ^ <(-,*) , 

2 - o o * 
a 

~ < Me 

~ AH'* + Mt + M 

o 

where AM' is of the form (42), except that it now includes an extra 

• • • * -2 

terra proportional to [ ( 2n^p+q) ( a>fM)~qaH*m0^] in the coefficient Mnmpq 

fV ■ 

(corresponding to s ^ nmp q in (43)), because of the double integration 
of Aa^®\ while AM^ is such that the inertia condition AM^®^(t^) = 0 
is fulfilled. The terra M > t can be included in the secular component of M 
(explained in part (b) of this paragraph) , so only the periodical part 


5 7 


“' (8><t) ' nmlipq(nonres) nma ^ nin P^ COS ^ n "'^ P+ ^^ 


+ < t > - •?) + am* V 

nmapqo 2 o' 


( 42 ") 


shall be considered here. The coefficients s £ nm pq in (42) and (42') are 
proportional to the Inverse of the frequency. Accordingly, the closer 
this frequency is to zero, the larger and, thus, the 

contribution AO § of AC to the As". . At about 1000 km height 

an orbiting spacecraft completes between 13 and 14 revolutions in one 
"nodal day 0 , which is the time - J 2 tt/ 0 " | it takes the. Earth to do a 
full turn with respect to the precessing orbital plane. If this number 
were exactly 13 or 14 then, calling it Nr and assuming that fj =* 0, there 

would be always some combination of n, p and q such 

• • 

that (n-2p+q)M + m0' * 0 for all potential coefficients of order 

m = k N , where k = 0, 1, 2, 3, . . . 

R 

/-w 

When the frequency is zero, s^^^ becomes infinite and the solution of 

the linearized equations given by (42) and (42' ) does not apply any 
longer. This case, to be discussed in the next paragraph, is known as a 

perfect resonance . In practice, is quite small for altimeter 

satellites, but not zero, If the number Nr for one "nodal day" is still 
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very dose Co an integer, the effect roust be very similar to that just 

described, at least over a period of time of a week or so, which is much 

smaller than the period of 6). The perturbations, though very large, become 

periodical, in agreement with (42) and (42') » even if their periods are 

very long. This is a case of deep res on ance , and the resulting 

perturbations, if given enough time to grow, may have peak amplitudes of 

the order of a kilometer or pjore; in that case, a linear theory is no 

longer adequate and one has to use a nonlinear approach. This is not the 

case for the errors of weekly altimeter satellite arcs, which never reach 

more than a few meters in size. In reality, Nr differs considerably from 

an integer, though those AC errors whose orders come closest to 

run a 

satisfying the resonance condition given above, or resonant orders , still 
can produce long period effects much larger than those of the other 

orders. This is known as shallow resonance , and altimeter satellites have 
orbits with such resonances at orders close to whole multiples of 13 or 

14, In particular, when k * 1,2,3 one has the so-called primary , 

secondary and tertiary shallow resonances, respectively. Orbits 

where H) - 0 and there is an integer number of turns Nr over a whole number 
of -nodal days" Nr, are known as "frozen", repeat orbits . The same 

argument used when Nr = 1 applies now: these orbits must experience 

perfect resonances with those AC nma that satisfy the condition m =* k Nr, 
which includes the zonal s for k = 0. This type of orbit turns out to be 
quite important in satellite altimetry, and more shall be said about it 
presently. 
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Next in importance to perfect, deep or shallow resonances are the 
effects »£ those terras in (42) and (42”) where (n-2p+q) * 0, as the 
frequency then is reduced to -qu + mV, For the orbits under 

consideration, this sum is quite small compared to the "once per 

# • • 

revolution" frequency, approximately equal to (M+u>) M, so it must be 
much smaller than any other where (n-2p+q) * 0, provided that q and m are 
small integers. At altitudes of 1000 km, |q| <2 for all terms of any 
real account, and the perturbations are substantial only if m < 10. Given 
the prevalence of mV over q&>, these frequencies and the corresponding 
terms are usually known as "m-dailies" . Because of the size of their 
effects, which show quite clearly in data from satellite tracking, deep 
resonances, shallow resonances and "m-dailies" are quite useful for 
estimating the AC nma from this type of data, as they comprise most of the 
signal. For the same reason they are quite important in understanding the 
nature of orbital errors, which is the subject of this work. 

The total amplitude of the oscillation associated with any particular 
frequency in (42) and (42”) is 

l AC s . , where (n-2p) , q and m must be all constant, 

L nma inmpq ’ v ’ ^ 

npq ■ 

so only AC of the same order m and degree n of the same parity (even, 

odd) can contribute. This amplitude, being a weighted sum of ^ nmu > 

known as a "lumped Coefficient" ; The weights in the sum are functions 

of inclination and of eccentricity through the F (I) and G (e) 

? nmp npq . 
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functions in expression (38); which AC dominates the sum will depend on 

nmot 

both T and e* Thus, modelling the gravitational field with tracking 

data requires satellites with a wide variety of inclinations and 

eccentricities, to disentangle better the various AC ; also for a given 

nmO 

order m, coefficient errors of e /an degree n produce perturbations most of 
whose frequencies differ only in ± qfi> from those of odd degree (the main 
ones having |q| < 2) , so it is desirable that the interval between the 
first and the last observation of a satellite be at least a substantial 
part of its apsidal period 2 tt/£>. 

Of all the C , the zonals are the best known at present; their 
nmot 

determination involved most of the early work done in the field of mapping 
gravitation with artificial satellites. It was the analysis of the 
tracking of one of the first USA orbiting spacecraft, a VANGUARD, that 
revealed the existence of a strong third zonal indicating that the Earth 
is rather "pear-shaped", with the southern hemisphere slightly larger than 
the northern one (Eckels et al. , 1959). 

(b) (n-2p+q) (M-ko) - qCfrhnB'* = Q (Resonance) 

As already explained, this case is known as a perfect resonance: it 

occurs whenever the angular frequencies M, w and 0' are locked in step 
with each other, which is a rare case, but quite relevant to the study of 
the orbits of altimeter satellites. Then, the orbit crosses the 

gravitational perturbations associated with some of the AC at 


61 


intervals equal to the period of M, the natural frequency of the 
linearized system that, forced in this way, resonates . One example is the 
"frozen", repeat orbit , whose ground track repeats with a period of N Q 
"nodal days" which is an exact multiple of that of M. An extreme case of 
this is the orbit of a geosynchronous satellite (often used to relay 
communications) whose ground track, ideally, is on a fixed point on the 
equator where it "repeats" for ever. Here the main resonant coefficients 
are C 220 » a nd C 221 * related to the "triaxiality", or lack of rotational 
symmetry, of the best fitting ellipsoid for the Earth, 

As given here, the resonant condition requires an excitation of 
0 frequency, instead of M. This apparent contradiction with the 
physical explanation given just now is resolved by noticing that whenever 
(n-2p+q) (M+o>) - qaH-m0'=O then (n-2p+q±l) (M+w) - (q±l) uH-m0 , '= ± M, so 

any ^ nma that produces a zero frequency term must produce also forcing 
terms of frequencies M and -M (either of which makes the linearized 
system resonate). Moreover, the amplitudes of the 0 and ± M forcing 
terras are implicit functions of each other, as both depend explicitly only 
on the parameter q (in all of them, a., 1 , 5 , n, m and p are the same). 

As before, the linearized equations for Aa and he can be solved 
independently by direct integration with respect to time. Adding up all 

< (or) 

resonant terras and symbolizing their sums by As^ , 



{ y AC s . cos ( $ +■?■)} (t-t^) 

* ^ nma inmpq v,f nmapqo 2 1 ft 

nraapq(res) 


< , , A <(g)k 

s.t + As. ■ 
i io 


( 44 ) 


where " I " indicates the sum over all the resonant combinations of 
nmopq(res) 

n,m,p,q, while Agj: g ' ~ “ ®i C o ( so inertia is not violated at the start 

of an arc), and ^ is the expression in curly brackets. The equation for 

AM, on the other hand, has a solution that depends on Aa because of 

the not negligible "interaction term" - — Aa: 

2 a 


{ I AC M cos(<j> „ )} (t-tS 

1 • L nma nmpq nmapqo o 

nmapq(res) 

“ | Z + Al^ g)k (t-t^)] . (45) 

a 

In practice, the quadratic term in (45) is unlikely to become large enough 
to matter within any weekly arc. However, quadratic terms have to be 
Included not only in AM, but in Aa and Ae as well, because true 
resonances are deep rather than perfect, with frequencies that are very 
small but never zero. The corresponding terms have, therefore, very long 
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periods, which over a much shorter interval, such as one week, can be 
approximated by quadrics: 


A <(g) *k „2 , <k , k 

4s 9 w s t + s t + s . 
i i i io 


Including both periodical and deep resonant terms (if there are any), 
the complete solutions of the three linearized equations of motion (38) 
corresponding to gravitational model errors can be written as 


ia (8) - ^ (s) + ?t 2 + $ k t + a k 


(46a) 


~(g) ‘k2 , <k k 

Ae +et + et + e 

o 


(46b) 


= iS (g> + SV + i* k t + M k 


(46c) 


k k. k 

where the constant terms a , e and M are such that, at the start of any 

o 9 o o 

weekly arc, Aa^^(t^) = Ae^(t^) = = 0. The second degree 

terms in (46) are likely to be considerably smaller than the linear ones, 

(p) 

so the non-periodic parts of the As ^ may depart only slightly from 
straight lines* 


Errors in drag, etc., can also have very long period effects 
resembling weekly quadrics, as well as others looking rather like higher 
degree polynomials (expression (18)) • 



*v 


3.14 Orbit Error Due to Incorrect GM 

At present, p (often referred to as "GM", where "G" is the 

universal constant of gravitation and "M" is the mass of the Earth) is 

known to better than six significant places. According to (26), the 

effect of an error Ap is equivalent to an error C ^ in each 

|j nm ot 

potential coefficient, or about C x 10 - ^. This means that the only 
v nmcx J 

appreciable effect is that of -jp c qqq* as the zero degree (or "central 

force") zonal (which has a unit value by definition) is almost three 

orders of magnitude larger than C 20 O* and five or more than all the 

rest. Therefore, to understand the influence of Ap on the computed 

orbit, it is sufficient to consider the case where the field consists of 

the zero harmonic alone. In this field the orbit obeys Kepler's laws, 

. ’ 3. 1/2 

so its frequency is M = (p/a ) . For a nearly circular orbit of 

mean radius r, a ~ r, so M ~ (p/r ) ' . An error Ap in p requires a 

compensating bias Ar^ in r to keep the frequencies of the true and the 
computed orbits very close; accordingly, Ar^ = • Otherwise, the 

along track errors may become very large after several days. While Ap, 
through its influence on the adjusted initial state, causes more than a 

bias in Ar (see expression (23)), here Ar^ is, nonetheless, its most 

distinctive and important effect. This bias is virtually the same for all 
the arcs, as long as r does not change much. Such is the case throughout 
most of the mission of an altimeter satellite, where the effect of drag is 
relatively small because of the height, so Ap will cause the 

alt imet rlcaily determined mean sea surface Tj , , for instance, to appear 
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%'■ ■ J*. 


higher or lower by the constant offset Ar^. Likewise, there may be an 
error in the size of the reference ellipsoid, approximately equal 
to AR (which may be partly related to Ap) , so the estimate of the mean sea 
surface topography w will be biased as well by Aw = r/3( Ap/p)-AR * 

°C O 

However, for any meaningful oceanographic application, what matters is the 
slope between the sea surface and the best fitting equipotential; the non- 
zero global mean value (bias) of the computed topography can be ignored. 

3.15 Expression of the Radial Error Excluding its Tidal Part 

Going back to (16), one can see that Aa appears in Ar directly, 
whila Aa, Ae and AM modulate together an oscillation of frequency M. 

This causes the frequencies present in (42) to appear in Ar both 

• • 

unchanged and shifted by ± M. Because of the -M shift, terms whose 
original frequencies are very close to M produce very long period 
oscillations in Ar, resembling, for a weekly arc, quadrics that can be 
lumped with that from the direct contribution of Aa. The +M shift, 
on the other hand, converts the sum of the deep resonant terms 
in Aa, Ae and AM into an oscillation of frequency $ modulated by a 
quadric. 

Putting together the various parts that make up the radial error, one 
can get an approximate "weekly" expression for Ar. To this effect, 

replace (46) in (16) while taking into account (18), (23), (42), (43) , 
(45) and the usual trigonometric identities for the sine and cosine of the 
sum and the product of two angles. The result is: 


66 


A , s< NT > V 

Ar(t) = ), 


• • • • 


nmotpq(nonres) 


AC {r cos [ ( (n-2p+q) ( — q tiHTn ©") ] 

nma nmpq nmotpqo 


+r cos [ ( (n-2p+q+l) ( wfM)-(q+l ) o>fm Q”) t+(j) ] 

nmp(q+l) nm«p(q+l)o 


+r cos[((n-2p+q-l)(uK-M)-(q-l)tiH-m*e')t+(j) ]} 

nmp(q-l) nmotp(q-l)o 


+A cos Mt+B sin Mt+t+C t cos Mt+D t sin Mt 
k k k k 


2 , 2 

+E t cos Mt + F t sin Mt 
k k 


+ j 0 • 


(47) 


where J is a small integer, t Q is such that M q = 0 in (40c) (to 

eliminate a term (q ± ^) M q in the arguments of the cosines and simplify 

the formulae) and the r n depend on the AC * Notice that the 

nmp(q ± '/) ™“ 

subscripts "(g)" and "(ng)" have been replaced by "(NT)" (for "non'" 
(NT) 

tidal"), as Ar N has here both gravitational and non-gravitational 
causes, but excludes the effect of tides (to be discussed later). The 
coefficients r nmp q are 


r = a 
nmpq nmpq 


r nmp(q+l) ~ 2 -nmpq 


= 7 r (M a e - a e-e a) 
o nmpq nmpq 


- 1 ) 


- tt (M ae + a e + e a) 
2 nmpq nmpq ntnpq 
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so they are functions of a, e and X (through the F nn)p (I), G np q(e), etc., 
in (38a-c)), but not of M, cu or Q. For a satellite like GEOS-3, where the 
mean ellipse changes appreciably over the mission due to uncontrolled 
orbital decay caused mostly by air drag, the r nmp q, r nrap (q+i) and 
6 , d> / .i\ are also arc dependent, although the former are not 

likely to vary appreciably from week to week, This decay was prevented by 
orbital corrections in SEASAT and shall be prevented likewise in most 
future altimeter satellites, so expression (47) applies best to them. For 
arcs of up to one month (47) is still likely to be valid if the envelope 
of the increasing ft oscillation is approximated by a cubic or a quartic 
instead of a parabola. 

The resonant part C^t cos fit + D^t sin fit + E^t^cos fit + F^t^sin fit 
caused by gravitational field errors is only important in deeply resonant 
orbits like those discussed in the next chapter where, as shown by 
computer simulations, i£ tlva AC were of the order of the published 
accuracies for the coefficients of GEM 9 (Lerch et al. , 1977), these 

increasing oscillations could build up at a rate of several meters per 

week. As for the periodical terms, the power in the spectrum of that 

part of Ar caused by errors like those in GEM 9 would be distributed as 
follows, in the ce-sg of SEASAT: up to 1 cycle/revolution, 3m (r.m.s. ); 

from 1 to 10 cy, t lev. , 2m; 10-20 cy./rev. , lm; 20-30 cy./rev. , 0.1m; 30- 
40 cy./rev. , 0,05m; above 40 cy./rev., less than 0.05ra. More recent 

gravity field models, such as PGS-S4 (Lerch et al., 1982), are thought to 
reduce considerably all the components of Ar associated with them. 
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4 . 0 flggKE S- JOTB RKPMTIMG ABOUND TRACKS 

4.1 Tfcf yytozen i,r . Repeat (Orbit 

For -wmf ^pplfcillcsls of altimetry it is desirable to choose the 
ohMit so sattellite passes over the same places on Earth every few 

diifS# repeating its groand track periodically. This was done, 
©spexfslftOtall'3^ during the lassfc 'month of operation of SEASAT (Cutting et 
rfX#* . 1978),, and 1 is gilanraed for all civilian altimeter missions now 

.J>eirSg (.considered. Because of the disruption of the chosen orbit by drag, 

" grs* : sure» and other disturbing forces, the ground track cannot 
te rep 'ated exactly, but it is sufficient to keep it within a band a few 
kilometers -wide (about ±5 km for SEASAT (Bauer, ib. 1978), and probably 
abont *1 km for the satellites that will follow over the next decade). 
This can be achieved by firing small rocket engines to correct the orbit 
every &o often, for example less than once per month for ERS-1 (Dow and 
Klinkrad, ib. , 1982) , so as to turn the drif t back towards the ideal 

orbit. The maneuvers must be brief, as the thrust of the rockets cannot 
be modeled accurately enough for precise orbit determination, so 
measurements taken while they are operating are hard to interpret. For 
this reason, the precisely estimated ephemeris cannot run through one of 
these maneuvers; an arc must end just before and another one begin just 
after it. To get the desired repetition, it is sufficient that the 
maneuvers keep the mean ellipse of the estimated orbit so that, if this 
ellipse were the actual trajectory of the spacecraft, the ground track 
would repeat perfectly. For this, such mean trajectory should start and 
end at the same point (in Earth-fixed coordinates) in each repeat cycle. 
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ERS-1: Projected 3-day repeat groundtrack 


ORJGWNAL' PA' V • 
Of. poor quauthc 



.AVAwAl^VAWAWAVAWAVAWAVAMAWAWAMWAUAUAVAVAWAWAUAWAWAUf, 

VivAftvXwAVAVAVAVAVAVAVAVAVAVAViffiVAVAVAVAVAVAVAVAVAV.AVAVAVA* 

^M^M^iVA^MAMAMAV.MAVAMAMAMAwAVAVAMAMAVAViMAMi^VA^wAvAMAVAMM 


mmm smmmmmm mmMmmmmmmmmmm m 

iiWiiiiW (Ni 

mililmU 

wmmmro jg |iA«o^^ 



mmmmmmi 


POSEIDON: Projected 26- ^ay repeat groundtrack 
(shown only over the North Atlantic) 


FIGURE 2(b). Repeating Groundtracks of ERS-1 and POSEIDON 


71 




This is possible only if w, ft and &' ■ ft - & are all exact multiples of 

the fundamental angular frequency w ■ 2 n/T , where T is the 

rc re rc 

repeat period. For altimeter missions T rc has been chosen, up to now, of 
the order of a few days (three days for SEASAT, nine for TOPEX, three for 
ERS-1; Figure 2 shows the ground tracks of these satellites). The 

approximate expression (11) for Jo would require a very large mean 
eccentricity of the orbit (close to e * 1) for the perigee to precess in a 
few days, making &T rc • k2 7t ( k»0 ,1,2...), except when is zero. 

Unfortunately, such an eccentric orbit would have its perigee deep 

underground, This leaves 5 j ■» 0 as the only practical choice. The 

question is: can one have both a small eccentricity and an Jo that is 

virtually 0? All zonals can produce secular perturbations in u>; formula 

(11) gives only the effect of the largest of the zonals, C^qq* However, 
when e is very small, the position of the argument of perigee can be 
rotated through a large angle by very slight changes in the shape of the 
orbit, like those produced by the other zortals, particularly the second 
largest, C 300* As these perturbations are of different signs, their 
combined effect on the mean ellipse may Cancel out that of C2 jq alone, 
thus making 5) = 0 regardless of the inclination. Cook (1966) developed 
a theory for the mean ellipse to take into account the effect of an 
arbitrary number N of zonals. From Cook's theory one can get the values 
of e - and Hi that "freeze" the mean orbit assuming that "a,T and the first 
N zonals are known: 
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e 


t N(odd) 

"* J 3 °"°° 

* 1.182 x 10“ 3 


(l)" ^ 0 Kl^ e T >1°200 (I)' 

(|) sin I + 0 (C 500 ) , 



(49) 


where the summation in (49) is over odd values of n. Notice that (48) 
fixes the mean perigee at the highest latitude reached by the ground 
track. This theory involves a number of approximations and breaks down 
for 1 close to the critical inclinations (~63.4° and ~U6.6°); in this 
special case, the theory of Hough (1981) shows the existence of "frozen" 
orbits with small e and heights of about 1000km where w = as in (98) , 

although the equilibrium value of e is not given by (49), but by another 
expression that also depends on the first N zona Is. TOPEX's orbit is 

likely to have an inclination near the critical value of 63.4°. 


Cook's formulas are used to this day (for example, see Dow and 
Klinkrad, ib. 1982). It is important to notice that expression (49) is a 
continuous function of "a and T, so whenever a "frozen” orbit is possible 
at a given height and inclination, there are infinitely many others at 
heights and inclinations close to that one. SEASAT is a clear indication 
that this is the case at typical heights for altimeter satellites. How 
well can one "freeze" an orbit in a low degree zonal field? Recently I 
integrated numerically "frozen" orbits at various inclinations, using as 
the only driving force that of the field of the 9 first zonals of the 
GEM 9 model (Lerch et al. , 1977). The initial conditions were adjusted 

iteratively, starting with the values given by (49), until, making 
allowances for the precession of the orbital plane, they "returned" to 
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their starting positions and velocities within 10“^ meters and 10*"^ meters 
per second, respectively (Colombo, ib. , 1984), Of course, a real orbit is 
subject to numerous perturbations that prevent it from being that "frozen" 
and, in any case, the accuracy with Which a satellite can be put in a 
particular orbit (and kept there) is also limited. Fortunately, the kind 
of accuracy that can be achieved is enough for altimeter missions. 
Moreover, as any large perturbations are also likely to be slow in growing 
(i.e,, have long periods), they can be kept in check by occasional orbit 
maneuvers, 

4.2 Periodicity and Resonance 

A complete first order perturbation theory for "frozen", repeat 
orbits is given in (Colombo, ib* , 1984). What follows explains the part 
of this theory relevant to radial orbit errors. 

When w => 0 and the satellite completes a whole number of 
revolutions in an exact number Nq of "nodal days", the frequencies of 
the trigonometric functions in (38 a-b) and in (47) are all of the form 

(n-2p+q) M + m0' = ((n-2p+q) N - m N ) w , (50) 

K u rc 


"nodal", day .(lepgth equal to 2 7T/ 1 0 1 ) is .very close to one sidereal 
day because & in &' - v is much smaller than 0 . 
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where w rc is the angular repeat frequency of the orbit. The numbers Nr 
and Np must be relative primes , i.e., have no integer factor in common 

other than 1. Otherwise, if the largest of such integers is N'g, there 
will be a repeat every days as well as every N D days, with N'g < Ng. 

The minus sign of ~mNg comes from 6’' * ft - b ~ -6 , so &"* < 0, while 

here &> rc is chosen positive, and M is always positive, according to (9) 
and (9'). As before, terms in the analytical expressions of the 

gravitational errors can be separated into two classes: non-resonant and 

resonant. 

The most important characteristic of the non-resonant part of the 
orbit error is that it is now a Fourier series whose terms are harmonics 
of the repeat frequency <o . Therefore, this part, which is the most 

complex component of the error caused by the field model, is a periodical 
function of time, and its fundamental period is the same as the repeat 
period of the orbit. 

This strong property reflects the rotational symmetries of the mean 
repeat orbit and its ground track, each coiled around the Earth in a 
spherical helix that closes on itself. 

Resonances in a ’’frozen", repeat orbit occur when the frequency in 

(50) is zero. This happens for all AC whose order satisfy 

nma 

m = k N r (k = 0, 1, 2, ,..) . 
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(51) 


The corresponding force errors are "encountered" by the computed 

satellite position at repeated intervals of exactly one revolution, which 

for a "frozen" orbit is also the period of ft, so the linearized equations 

of motion are driven by these errors at their natural frequency ft and have 

resonant solutions. For repeat periods of three days or longer (those of 

all satellites considered here), the lowest degree n of a AC that 

nma 

causes resonance must be n > in = k N > 43 (= 43 for SEASAT, with a 3-day 
period )* Expressions (32) and (33) show that , as n increases, the 

strength of the harmonics in. the expansion of V must eventually fall 
off. At 1000 km, even resonant coefficient errors with n ~ 43 should have 
a very small effect on Ar. So the only case where resonance is likely to 
be important is when k = 0 and, thus, m - 0 (zonals).C^ Expression (38a) 
shows that Aa cannot have a secular part produced by zonal errors, 
because Aa is proportional to (n-2p+q) and must be zero at resonance 
if u*=m=0. Therefore, field model errors may cause significant secula.. 
departures from the true values of all the elements but the semi-major 
axis, in contrast with surface model errors that can change this element 
as well . 


The main errors are due to the odd zonals, mostly because the term 
proportional to 3G np q(e)/3e in (38c) is very large when e is small 
and q =1 (3G npq (e)/3e^2 1) , and causes the affect of AM to prevail 
in Ar over those 1 of Aa and Ae. If q-1, the resonant condition for m-0, 
i.e., (n-2p+q)=0, requires n to be odd. 
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4.3 The Stability of a ’’Frozen 11 Orbit 


In a "frozen" orbit, zonal resonances will occur regardless of 
whether the ground track repeats itself or not, because the condition for 
resonance (n-2p+q)=0 does not involve M or fr'j zonal resonance is 
inherent in any "frozen" orbit. Zonal force errors are not affected by 
the rotation of the Earth, which shapes the ground track. Since the known 
values of a, I and the zonal coefficients used in formula (49) may be 
accurate but not perfect, and the effects of higher degree and order 
potential coefficients, as well as those of non-gravitational forces (drag 
in particular) are not considered in (48) and (49), the supposedly 
fixed e and to will, in fact, change very gradually from their starting 
values. 

At any time, the mean orbit defines a point in a six-dimensional 
Euclidean space with Cartesian coordinates a, e, I, to, fi and M. As the 
orbit evolves, the point follows a slow trajectory in this space. Let Q 
be the point's projection on the plane of the axes e and to. If Q is 
initially identical with the "frozen point" P defined by (48) and (49), it 
will remain there indefinitely. Otherwise, as is shown in Figure 3, 
Cook's theory predicts that it will follow either a closed trajectory 
around P (libration) or else an open one where w eventually ranges from 
0° to 360° (precession). Similar predictions are made by Hough's theory 
(ib, 1981) for near circular orbits close to the critical inclination. 
The trajectory will be closed if Q is initially near enough to P, open 
(otherwise; in either case (taking Hi modulus 360°) Q returns to the same 
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starting position after the same length of time , which is the period 
corresponding to ui as given by (11). Therefore, the closer Q is to P, 
the smaller Its tangential velocity along its own trajectory, and the 
smaller, in general, the projection of this velocity on the w axis, or, 
approximately, the weekly average fi> (for SEAS AT Zo was always very 
close to 1° per year). 

Since the values of the zonals used to calculate the coordinates of P 

with (48) and (49) are only those of the N first C n QQ and are subject to 

errors AC ' , the true "frozen point" is P" and not P. As (48) does not 
nUU 

depend on any parameter, both P and P' are on the same line, parallel to 

the e axis, where 115= tt/2. Even if Q is placed exactly at P, it cannot 

stay there, being forced by the zonal errors AC to move instead along a 

nOO 

trajectory that encircles P', as shown in Figure 3. Initially, this 

trajectory will be very close to the local tangent, which is a straight 
line, so the variation in both e and "3 will appear to be secular when 
observed over a sufficiently short time. As P and P' are likely to be 

quite close, the movement of Q can be very slow and one week can be 

regarded as a "short time". Because of drag, the semimajor axis will 
decrease gradually, so P will Ishift upwards along the line W » ~ (see 
(49)), and Q will follow a slowly widening spiral instead of a closed 

curve. 

In the event of the mean ellipse becoming too different from the 
chosen one, the rockets in the spacecraft can be used to correct the 

orbit, bringing it back to the right altitude and also shifting Q 
to Q" (Figure 3) , on a path leading back (approximately) to P. 
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e p . - 1.182 x 10 — J (R/a) sin I 


FIGURE 3. 


Long period evolution of the "frozen" orbit in the e.oTplane. P' is the true 
"frozen" point, P the computed "frozen" point; the line from Q to Q' 
indicates a corrective maneuver after Q has moved too far away from P. 





4.4 T he Repeating Ground Track 


As clearly shown in Figures 2a-b, the ground track of a "frozen", 
repeat orbit is a curve of strong symmetries, including: 

1. rotational symmetry about the z axis; 

2. periodicity (or self-closing, ot time symmetry); 

3. equatorial symmetry of its crossover points (approximate, the 
same as (5) below, because of the slight eccentricity of the 
orbit) ; 

4. mirror symmetry about each of a set of 2N^ equally spaced 
meridians; 

5. equatorial antisymmetry of the values of the At (the intervals 
between ascending and descending passes at the crossovers), 

(a) Geometry 


As a consequence of (1), (3) and (4), the crossovers are arranged 

along parallels; there are N R equally spaced crossovers along each 
parallel; these parallels come in pairs equidistant from the equator 
(which is one of them only if N R and Njy are both odd numbers). All 
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crossovers are aligned along the 2N R meridians of symmetry and their 
connecting arcs form diamond-shaped partitions, each with the left and 
right vertices on two alternate meridians and the north and south vertices 
on the intermediate one. 

Regarding the points of extreme north and south latitudes as 

"crossovers" where At=0» there are at most Nr .+ 2 parallels with 

crossovers. Over the northernmost, perigee is reached, apogee over the 

southernmost; there ascending passes end and descending ones'* D begin, 

or vice versa. At each parallel, the intervals At (Modulus T rc ) have 

two different values of opposite signs: At and At, , where 

a b 

| At a | + ( At ^ | = T rc ; this is a consequence of (2). If Nr and Np are 
both odd, the equatorial crossovers have | At a | = J | = N R T Q /2 = NpTp/2 

= T rc /2 (where T Q is the orbital period, and Tp the length of the nodal 

day). If either N r or Np is even (both cannot be even simultaneously, as 
they must be relative primes) the equator cuts through the middle of the 
lower half cf one diamond, the top half of the next, and so on. If Np is 

even and Nr is odd, there are at most Nr + 1 parallels with crossovers. 

The mirror symmetries about meridians with crossovers can be understood by 
imagining that, as the satellite reaches a crossover along an ascending 
pass, the North-South component of its velocity vector is suddenly 

reversed: In Earth-fixed coordinates the satellite (and its ground track) 

will follow a course symmetrical to the one that brought it to that 

crossover; this new course (when m - ™) is that of the descending pass 

(Dfj < 0 and R > 0, respectively. 
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over the crossover. The equal spacing In longitude between meridians with 
crossovers is a consequence of both this mirror symmetry and of the 
overall rotational symmetry. The equatorial antisymmetry of the 
intervals At follows from the approximate antisymmetry (respect to the 
geocenter) of the Earth and the slightly eccentric orbit when they are 
orientated according to their respective senses of rotation and 
revolution. 

From the figures it is clear that the crossovers tend to crowd 
together towards the northern and southern edges of the band covered by 
the ground track. The overall pattern is that of a fishing net with its 
stitches aligned vertically and stretched in that direction towards the 
middle, without horizontal shrinking. 

(b) Finding the Crossovers 

For a crossover to take place, both the Earth and the satellite must 
come together twice: first during an ascending (or descending) pass, and 

again during a descending (or ascending) one. This means that, at the end 


of an interval equal 

to At, both the 

latitude 

and 

longitude of 

the 

subsatellite point must 

be the same as 

they were 

at 

the beginning. 

It 


also means that, since the argument of perigee is "frozen" at to = — , 
the orbit must follow a path symmetrical to the meridian plant of the 
perigee in a system of coordinates where this plane is fixed. The 
beginning and end of this path occur at the same latitude $ both in this 
system and in the Earf.h~fixed one, and the time it takes the satellite to 
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complete it is At. Therefore At is a function of <j>j it also depends on T 0 
and on the orbital inclination T, so At *= f (<f>, 7, T^). This relationship 
and the two ground track equations for <j> and X (Table lc) provide a set of 
three equations in the unknowns At, <p and A. Unfortunately, these are 

transcendental equations and have, as far as I know, no closed-form 
solution. Without this, to find the crossovers and their At it is 
necessary to search for approximate solutions, exploiting the symmetries 
of the ground track to save effort. 

Starting at the highest latitude <{>=7, choose a point in the ground 
track (!}) o ,A o ) directly under the perigee (i.e. , ^ q = 7) ,so the corresponding 
true anomaly is M, Assume that the time of passage of the satellite 
over this point was t Q =0, and shift the origin of longitude so A q =0 , which 

— Tf 

makes 0 = fi = - — (in this paragraph A is defined, for convenience, in the 
interval -it<A<tt ) 5 select a small time-step T « T Q . The line A = A q is 
one of the 2N R meridians with crossovers. To find the others, use the 
equation of the center to calculate f^. at t=T (see Table lc) knowing that 
M t = 11 T, where n is the mean motion, and obtain <j> T and A^ with the ground 
track equations. Repeat this process at t=2T, 3T , . . . , kT, continuing 
along the same descending pass until reaching a point ( , A^) just past 
the first meridian where | A| > ir/N R . Estimate the location of this 
meridian and the time tj of passage through it by assuming that the 
satellite was moving with uniform velocity relative to the ground between 
points k-1 and k. If necessary, refine this estimate by starting from 
point k-1 with a time-step smaller than T. Once the location of the 
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intersection with the meridian is sufficiently well established, the 

coordinates of this intersection correspond to a crossover point, 

and the other crossovers along the same parallel have coordinates 

(<J> 1 , + j Tr/ty, j»l,2,,..N R _p Find next the coordinates of the 

ground track points at times tj+T, tj + 2T, . , . ,ti+kT, until j >.| > 2it/N R , and 

repeat the search procedure just described, to locate the meridian where 

| | 2tf/N R and the corresponding tg, y, \ ° * a seconc * crossover. 

This determines, once more, the positions of all crossovers of 

latitude $ , Continuing in this way, it is clear how the remaining 

crossovers (<j>. , K.) and their corresponding times t., are to be found, 
i l x 

When this is done, the positions of all the other crossovers in the ground 
track are also known. To obtain the crossover intervals between the 
ascending and the descending passes, for each value of i take At' = 2t i 
if the satellite motion is prograde (same sense as the Earth's rotation) 
and At' = -2t if it is retrograde . Next, find N R ^, the smallest positive 
integer for which 


2|\l = (2 tt N q (N Ri /N R )} (Modulus 27r), 


The values of the two crossover intervals common to all crossovers at 
latitude ;<J^ are | At fli | = N Ri T Q + At' ( At gi > 0, if the motion is prograde. 


At a ^< 0, if retrograde), and 


At 


= T - At . I , the sign of At, 


being opposite to that of At 


ai 


bi 1 rc 1 ai 1 bi 

When great accuracy is not needed, one 
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may find the crossovers in the northern hemisphere, down to the equator, 
and then take advantage of the near equatorial symmetry of the ground 
track to locate the southern ones approximately, as the mirror images of 
their northern counterparts, using then the expression 

t (N R ~i+i) “ V 2 ~ H 

to finish the procedure in half the time. 

Finally, to put the crossovers in the Greenwich longitude coordinate 
system, add the actual geographical longitude of the starting 
point, X Q , to each of the A^. 


4.5 Orbit Error in a ’’Frozen 1 ', Repeat Orbit; Its Geographical 
Characteristics 


Neglecting, for the time being, the gravitational effect of the 
tides, the expression for Ar of gravitational and non-gravltational origin 
is the same as (47), except for the frequencies in the cosine terms: 


■ V AC {r cos[((n-2p+q)N 13 -mN r .)w t + 4> w ] 
(t) u /nma A nmpq N 1 ■ R D' rc nmapqo' 

nmapq(nonres) - 


+ r , . , N cos [ ( ( n- 2p+q+l ) N -m N )u> t + <f> , n ] 

nmp(q+l) R D rc nmap(q+l)o' 


+ r n»p(q-l>“ sl( <"- 2 ' ,+,1 - 1) V" W* + + nm<.p(q-l)o n 


# 2 , 2 . 

+ A cosMt+B sinMt+C tcosMt+D tsinMt+E t cosMt+F t sinMt 
K k k k k k 


'+ I r ik 

j-Q. JK 0 


(52) 


Here, in the 4> fo , , the term '4 ( l-(-l) n " m )" of (39) 

nmap(q ± l*}) o _ 

can be replaced with — ( l-(-l) m ) , because now = — (Colombo, ib. , 
1984, Ch. 2, par. 2.2). The very small u> due to the libration of the 
orbit about the ’’frozen point” has been ignored, in the arguments 
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C(n-2p+(q±°))N R -mN D )w rc t, because its effect on the true value of these 
arguments must remain negligible over weeks and even months, as explained 
later. 

Since the mean orbit is "frozen", and kept so by corrective maneuvers 

wherever necessary, the r , , , which depend only on a, e and I 

nmp(qldb JJ) o 

among the elements, are virtually constant. So it is now both possible 
and convenient, as shall be seen when discussing altimetry in connection 
with ocean tides, to choose a time origin for the phase angles 
^nmap(q + j±j}) o ( ex P ress *° n (39)) that is different from that of the 
start of the arc t£ (which still appears in the arc-dependent part at the 
end of (52), either implicitly or explicitly). Instead, in the periodic 
terms of a "frozen", repeat orbit , "t 0 " is the first time when M « 0 at 
the beginning of the repeat part of the mission. 

The periodical part of the orbit errors, according to (52), must look 
"frozen" to an observer who occupies the same geographical position every 
time the satellite passes by, as both the errors and the sightings occur 
with the same fundamental repeat period. The linearized equations of 
motion (including those for I, ai and ft not given in (38) ) show that 
periodical parts are present in all As^, and therefore in all 
combinations of them besides Ar, like the along and the across track 
errors. So the differences (or residuals) between tracking data and the 
corresponding computed values of the ranges or range-rates (computed 
according to the positions and velocities of the estimated orbits) may 
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exhibit also periodical parte as they are, to first order, linear 
combinations of the orbit errors. 

After subtracting altimetrlc heights determined at about the same 
locations along overlapping passes, temporal variations in the sea surface 
must be observable unless they are synchronous with the repeat period of 
the orbit, but the periodical part of Ar must vanish. This is, of course, 
very convenient when one wishes to study such changes in height as tides, 
etc,, because a substantial part of the orbit error can be "filtered-out" 
in a simple way. 

Because of the existence of a large periodical component that depends 

only m position along the ground track, Ar must show a strong spatial 

correlation. This correlation was found, empirically, by Anderle and 

Hoskin (1977) when they analyzed computer-simulated values of Ar. 

Inspection of (52) also shows that, if Ar^ and Ar^ are the values of the 

periodical component along ascending and descending passes, respectively, 

then Ar and Ar\ depend on <j> and A alone. Moreover, Ar and Ar , are 
a d a d 

different functions of position, so their differences at the crossing 
points of ascending and descending passes are not zero (except for some 
components of zonal origin, as explained later). In the more general case 
(expression (47)) of satellites like GEOS-3 , a smaller part of Ar 
(corresponding to terms where q=m*Q, for example) has these periodical 
characteristics exactly, but still a considerable degree of geographical 
correlation can be expected for errors within passes separated in time by 
less than 1/4 of the apsidal period 2 tt/ to (i.e. , less than a month apart ). 
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How "frozen" an orbit must be for the effective cancellation of 

the periodical part of Ar in the differences between overlapping 

passes to take place? All significant components of the orbit error 

due to the field model, at a height of about 1000 km, are produced 

by AC with n,m < 40, so they must have wavelengths V longer than 

nma 

\i. 1T „ ~ 360°/40 <• 0°. The time T between orbit corrections should be 
MIN 

short enough to prevent an excessive departure from exact periodicity when 
w * 0. The unwanted phase-shift due to to, according to (47), is 

qU)T, where |q| j< 2 for any terms big enough to matter. Therefore ^ 
qo)T must be less than I v min (because a sinewave may reach its full 
amplitude in a quarter-cycle), so 

q° i 

T < f . (53) 

* 28 ) 

In the case of SEASAT, 8) ~ 1° per month, so T < 1.13 months. Orbit decay, 

mostly produced by air drag, changes a and, consequently, M, and 

through h (see expressions (9'), (10) and (40a)). This causes a growing 

misclosure of the ground track that can reach several kilometers within T 

months. This is much less than v , and the effects of variations in the 

MIN 

other mean elements are smaller still, so the motion of the perigee is the 
main factor affecting the repetition of Ar. While T is not going to be 
the same for all satellites, it is likely that a maneuver a month, at 
most, will be enough to keep the orbit satisfactorily "frozen" in all 
cases. 


4.6 Observability of Zonal, Initial State and Other Errors In 
Differences of Altimetrlc Heights 

As the time origin has been chosen when the "frozen" mean orbit 

reaches perigee, so M n = 0 and M = Mt, then <#> =0 in (39) for any 

u nmotpqo 

periodical terms in (38) which are associated with a zonal coefficient 

error (m = a = 0). Accordingly, expressions (47) and (50) show that, for 

the zonals , the periodical radial errors in the "frozen" orbit are sums of 

cosines of j Mt = jM (j = 1 , 2, 3, . . . ) , so they are even functions 

of M . With the perigee fixed at u = | , M is uniquely related to the 

latitude of the subsatellite point; therefore, the periodical components 

of Ar due to the zonals depend on latitude alone . This conclusion is not 

affected by the choice of time origin, as it might appear at first, 

because the periodical terms under consideration belong to the "steady 

state" response of the linearized equations, which does not depend on 

either the time origin or on the initial conditions. It is easy to show 

that this latitude-dependence is true of all periodical errors associated 

with AC that are even functions of M and whose order satisfy the 

nma 

resonant condition (51), although those coming from the zonals are by far 
the largest, as already explained. 

A similar conclusion applies to the cos fit term (in 47), so this 
component of the "once per revolution" part of Ar, due to the initial 
state errors (expression (23)), must be also a function of latitude only, 
at least within the same arc. Therefore, in height differences along 
overlapping passes and in crossover differences, observability is nil for 
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periodical errors of zonal origin, and very restricted for any other 
errors whose periodical components have frequencies that are harmonics 
of M (Wagner, ib, 1984), 

The observability problems considered so far affect only those 
components of Ar in (52) of frequency kM (k=0, 1 ,2 , . . . ) and where 

{af[l-(-l) n m ] }tt/2=0 (i.e., sine terms are observable, cosine ones are 
not). For total lack of observability in crossover differences, the 
necessary and sufficient condition is that the intervals At between 
ascending and descending passes at all crossovers must be exact multiples 


of 

some "basic" interval 

V: 

At . (where At is the 

min min 

shortest 

of 

the 

At): if At = At ■ , then 

’ B min 

all 

At are exact multiples of 

At , . 
min 

The 


unobservable components of Ar, if the general condition were satisfied, 
would be w & 2Tr/Atg and its harmonics ka) u (where k=0 , 1 ,2 , . . . ). Assuming 


that there is at least one repeat ground track where the condition is met, 

and that for it Nr and Np are both odd numbers, this grid must have 

crossovers at the equator with At = N_T n /2, as explained in a previous 

paragraph, where T D is the length of the nodal day. This Atg must be a 

multiple of Atg, the same as At^^. If the inclination of the orbit is 

increased, and the seraimajor axis adjusted slightly to maintain a 

repeating ground track with the same Nr and N D , then Atg must remain 

constant, in nodal days, but At m ^ n must increase. Otherwise, as can be 

shown quite easily from purely geometrical considerations, At . would not 

min 

be long enough for the satellite to pass over the same point of the 
Earth's surface twice within one orbital period. So the new value 
of At m: £ n > measured in nodal days, must be longer than before, while 


Atg stays the same: the congruence needed to have totally unobservable 

frequencies no longer holds. In consequence, total unobservability is not 

an invariant property of the repeating ground track. This does not 

preclude the existence of tracks where the At are close to meeting the 

necessary and sufficient condition. In fact, as reported by Douglas 

et al. (1984), that of SEASAT came close to having At that were all 

multiples of At^^. Because Nr and Np are both odd in this case and the 

satellite is retrograde (it moves contrary to the rotation of the Earth) 

one can show, on the basis of the geometrical properties of its ground 

track, that if all At, including Atg, were exact multiples of At , then 

At min = Nr/(N r +1) T q , where T Q is the orbital period. As N R =43 for 

SEASAT, At . = 43/44 T « 0.977 T , so a = 2 it/ At . - 1.023 M, which is 

min o o’ u min 

the value of given by Douglas et al. (lb). In reality, each At missed 

being an exact multiple of At by a small margin, according to the 

rain 

diagrams in their paper. In any case, all periodical components of Ar 
associated with the AC have frequencies (in the case of SEASAT) rather 
different from or its harmonics, as can be verified by making N R = 43 
and Nq= 3 in (52). The non-periodical part of Ar (same expression) has 
most of its spectral power concentrated within two narrow bands, one 

containing u)=0, which is always unobservable, and the other centered at 

* • 

the partly unobservable M and, thus, also near « = 1.023 M. 

Nevertheless, as simulation studies indicate very clearly, there is no 
problem in observing part of this component as a gradual, cumulative 
change in the size of the crossover differences over a succession of 
repeats of the ground track. On the other hand, if one were to treat the 
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lines within these bends ais independent from each other and to estimate 
their amplitudes and phases from those same differences by least squares, 
it is highly likely that the solution may turn out to be unstable. The 
way around this difficulty is to acknowledge that the lines are not 

independent, but strongly related to each other through the few arc 

parameters A^, B^, etc. in (52), on which they all depend. In other 

words, this part of the error should be modelled, not as a sum of a large 
number of steady oscillations, but as the sum of the few aperiodic terms 
of (52), Instead. The bias and the oscillation cos Mt, both 

discontinuous from arc to arc, are not fully observable in crossover 
differences. Their changes in amplitude can be detected, if ascending and 
descending passes belong to different arcs, but not their individual 
values. 

There may be further limitations to the observability of Ar in 

crossover differences because work with these is restricted, in essence, 
to oceanic crossovers. Over land, the roughness of the terrain makes the 
interpolation of heights between consecutive measurements unreliable, and 
crossover differences are always made out of interpolated heights, as it 
is most unlikely that two measurements would be taken exactly on the same 
spot. 


| 4.7 Numerical Simulations of the Error in a "Frozen", Repeat Orbit 

| Table 2 presents the results of a computer sensitivity study of the 

effects of zonal resonances on the envelope of the "once per revolution" 

oscillation in Ar. Shown are the values of: the "AC ", which are 

nma 
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really the differences between zonal potential coefficients of two 

gravitational field models (PGS-S3 and PGS-S4, Lerch et al. , ib. (1982)) 

from n a 2 to n=9; the published standard deviation for the errors of the 

same zonals in a well known model (GEM 10, Lerch et al., ib. 1977); the 

contribution of each zonal to the increase of the "once per revolution" 

oscillation expressed as an average weekly rate; the sum of these rates, 

or total rate. The results correspond to the "frozen" » repeat orbit of 

SEASAT over an interval equal to its repeat period, or 3 days. They were 

obtained by first integrating numerically the orbit with one of the two 

field models, and then solving (also numerically) the exact (i.e, coupled 

and time varying) variational equations (21b) to get the individual radial 

effects of the as accurately as possible. As expected, the 

contributions of the odd zonals to the errors were the largest (see 

footnote after equation (51)). Since the size of the actual AC should 
. nma 

be much the same up to degree 29 (if one believes the published standard 
deviations for the coefficients of GEM 10 to be both true and 
representative of the "state of the art"), the overall effect of zonal 
errors may be several times larger than the total given in the table. As 
the pelagic tides are supposed to be only a few decimeters in amplitude, 
this error cannot be ignored when studying such oceanographic phenomena. 
Regarding the secular changes in a, e and M (not shown in the table) the 
calculations indicated that A& was, by far, the one that influenced Ar 
the most. 
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TABLE 2 


AC 


Increasing radial oscillations caused by gravitational zonal errors 
nOO* 


Degree n 

oC n00 (GEM9) 
X 10~ 9 

AC (PGS-S3/PGS-S4) 
niJU 

x 10“ 9 

Rate in 
meters/weel: 

2 

1. 

-0.51 

-0.02 

3 

1 . 

-1.65 

-1.40 

4 

1 . 

1.00 

=0.02 

5 

2. 

0.50 

-0.17 

6 

2. 

-2.68 

0.05 

7 

2. 

1.16 

-1.84 

8 

2. 

-3.80 

-0.06 

9 

2. 

-2.65 

5.41 


Total Rate; 


1.95 


These results correspond to a sensitivity study of SEASAT's repeat orbit 
where the zonal errors have been taken as equal to the differences between 
the zonal coefficients of field models PGS-S3 and PGS-S4. The published 
standard deviations of GEM 9 for the same coefficients are given for 
comparison. The coefficients are fully normalized and dimensionless. 


A fully normalized coefficient is: C 


(n+m) I 


1/2 


nma 


(2-6 m0 )(2n+l)(n-ra)! 


nma 
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The theory has been tested also by full simulations of Ar, done by 
F.J. Lerch and his assistant, N. Weiss, at Goddard Space Flight Center, 
whose help is greatly appreciated. "True" orbits, based on the ephemeris 
of SEASAT, were integrated numerically using the program GEODYN (Martin et 
al. , ib. , 1976) with a given field model (the "true field") and used to 
create "range data" between "terrestrial stations" (whose coordinates were 
specified) and the "spacecraft". With a field model (the "incorrect 
model") different from the "true" one used before, "computed" orbits were 
fitted to the "data" by adjusting their initial states. The "radial 
errors" were the differences between the "computed" and the "true" radial 
positions. These simulated errors, over a period of six days, exhibited 
the characteristics that one would expect from expression (52): the 
presence of resonant terms with linear and quadratic parts, the former 
exceeding the latter by an order of magnitude, and shorter period 
perturbations repeating themselves with the same frequency as the ground 
track, every three days. For these perturbations, the departure from 
perfect repetition did not exceed 1 cm in any of the cases studied. Carl 
A, Wagner at NOAA has obtained similar results (personal communication). 
The repetition of the shorter period components was checked by studying 
the differences between errors at passes of the "spacecraft" over the same 
points along the ground track at 3 day intervals. The results suggest 
that expression (52) represents quite faithfully the main qualitative 
characteristics of Ar due to the imperfections of the field model 
alone. Figure 4a shows what a typical plot of the resonant part looks 
like over six days, or two consecutive repeat cycles of SEASAT. The error 
caused by the field model has. a component at frequency M that is cancelled 
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FIGURE 4. Simulated radial differences between "true" and "computed" orbits 
(with adjusted initial state). 

(a) Tyical "bow-tie" pattern produced by GEM 9-like zonal errors. 
SEASAT, 6— day arc (two 3— day repeats). 

(b) Error caused by the drag model (Jacchia 1965— Jacchia 1971). 


(c) Error caused by wrong r eflectivity coefficient (set equal to zero) 

in the presence of radiation pressure. ERS— 1, one full 3— day repeat, 
heliosynchrorious arc (from Wakker et al. (1983)). 



by an oscillation of the same frequency and opposite sign produced by the 
error in the adjusted initial state which, in this way (see expression 
(23)), "absorbs" part of the gravitational error. This gives the envelope 
of the resonance its "bow-tie 8 * shape, a name suggested to me by Richard 
Eanes . 

The distribution of air density at satellite altitudes has two main 
large scale components! a zonal one (atmospheric oblateness) and a bulge 
that follows the Sun and thus turns once per year around the Earth's axis 
in inertial space. The atmospheric models used to calculate drag 
represent mainly these characteristics. Since they do so imperfectly, 
their errors have also a zonal part and a residual bulge that turns with a 
very slow angular velocity. The levels of the heating and ionizing solar 
radiation and of the magnetic flux do influence the air density} 
therefore, the density models are in error to the extent that they do not 
account properly for these and other variables. If these factors remained 
constant, both the large scale density and the model errors would be 
invariant in inertial space, except for the slow rotation of the bulge, so 
the "frozen" orbit would cross them with a frequency very close to M and 
the density errors would set off resonances. In reality, both solar and 
magnetic activities fluctuate, modulating the density. However, if one 
followed the practice of adjusting the drag coefficient once every day 
when estimating the orbit, the effect of all these fluctuations could be 
greatly reduced, as suggested by one of the simulations of Ar done with a 
"true" atmospheric density and an "incorrect model" (the models in Jacchia 
(1965) and Jacchia (1971), respectively), both used with the actual 


recorded changeo in solar and magnetic activities that took place in the 
second week of September of 1978. The gravitational field was assumed to 
be perfectly known, and the true drag coefficient of the satellite to be 
constant, like that of a "cannonball". Figure 4b shows the sort of radial 
error that was observed. In it there is also an oscillation of frequency 
M with an envelope somewhat resembling a "bow-tie", superimposed to which 
there is a fluctuation of much longer period, that can be represented by a 
low degree polynomial. While the envelope of the oscillation is more 
complex than in Figure 2a, the "secular" trend still dominates. Without a 
daily adjustment of the drag coefficient, the "polynomial" part was much 
more pronounced. Therefore, expression (52) describes quite well the 
radial errors observed in all the simulations (as long as the drag 
coefficient was adjusted once a day when non-gravitational forces were 
present). 

Like the distribution of atmospheric density in the case of drag, the 
main sources of radiation pressure (direct sunlight, albedo reflection and 
re-radiation in the infrared from the Earth) are "static" in inertial 
space, following the slow, apparent annual motion of the Sun. Therefore, 
one would expect part of the error in a model of this force to "set off 
resonances" in Ar, much as the errors in atmospheric drag do. For a 
heliosynchronous satellite in a repeat orbit, the resonances must be most 
pronounced, because the angle between the orbit plane and the direction to 
the Sun is "frozen". This effect seems to be present in Figure 4c, taken 
from Wakker et al. (1983) , Which shows the differences over a 3-day repeat 
of ERS-1 between two simulated orbits, one with radiation pressure applied 
and the other without it. 
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5.0 TIDES 


5.1 The Tidal Forces 

The attraction at a point x on the Earth's surface by a disturbing 
body 1 * (Sun, Moon, etc.) whose center of mass is at another point jr* can be 
discussed in terms of the gravitational potential of this body at x 


V*(x) 


y* 


r*-x 




(cos <|£) 


(54) 


where # is the geocentric zenith distance of the disturbing body at x, r* 
is the distance between the body and the geocenter, y* is the mass of the 
body times the universal constant of gravitation G, and R is the mean 
radius of the Earth. Here our planet is considered to be a sphere and the 
disturbing body, because of its great distance, a point-like object with 
all its matter concentrated at the center of mass. These approximations 
are sufficient for the present discussion. The ellipticity of the Moon 
may also have to be taken into account when estimating the orbits of 
artificial satellites over very long arcs (several months). 
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As the system of coordinates Is geocentric, it falls along with the 


Earth towards the attracting body with an acceleration 

E* « Total attraction of body on the Earth / Mass of the Earth (M ) 

—a e 

« G M* (r*) 

( 55 ) 

which is also an approximation, as the Earth is a sizable body of nearly 
ellipsoidal shape and not a particle, so its center of mass is not quite 

the same as its center of gravity in the field of the disturbing body 

_ . • . A 

(Wahr, 1979, Ch. 2 ). For an Earth-fixed observer, E cancels out the 

—•a 

atfrart'i.on of the disturbing body at the geocenter. Away from this point, 

the observer notices an increasing force. In an inertial coordinates' 

system, the explanation for this is that the observer only perceives the 

difference between TV* and the acceleration of his own quasi-inertial 

system E*. To use Newton's equations of motion in the observer's system 
—a 

in a simple way, -_E* is treated as the gradient of a fictitious 
gravitational potential 

^ 2 
U* rR , /R v \ 

r tee + < vr ) cos *5' 

■-f l & P „0 <=»= • < 56 > 
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V*(x) - f (- E* ) . dr - - 

geocenter 


where . dj:" indicates the line-integral between points A and B of 

the scalar product of and the oriented line element dr. By doing this 

with the force exerted by each individual body and then adding the 

respective W*'s in the right hand side of the equations of motion (19), 

these become formally like the equations of an inertial system and can be 

integrated accordingly. This is often done when computing orbits in 

"inertial space" (as the saying goes, though this is not to be taken 

literally). Expression (56) is the line Integral of -E* along a path 

■ —a . 

from the geocenter to with p*/r* as the constant of integration (a 
constant in space , not in time) . Adding (54) and (56) gives the 

expression of V*, the tidal potential as observed in the geocentric 
system, 


V*(x) = V*(x) + V*(x) 
i — — a — 



P n0 < cos V 


(57) 


According to the addition theorem for spherical harmonics 


P n0< « K>~ 2 { (,+ U> 

m=0 


-1 (n-m) ! 
(n+m ) ! 


P (sin <{>)P (sin <|>*)cos m (X*-A) 
nm nm 


■ 2 ! (Sol P nm < sl " « X *- X > - 


m=0 


nmO 


(58) 
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where (4>, A) are the geographical latitude and longitude of _x, while 

(<{>*, A*) are those of the point on Earth from where the disturbing body 

can be seen at the zenith (so 4>* is the declination and A* the right 

ascension minus 9 of the body). Replacing (-^r) n+ * P (cos #) with its 

K no x 

r 

expression according to (31) in terms of the osculating Keplerian 
elements a*, e* f I*, w*, M* of the disturbing body in the quasi inertial 
equatorial system , while truncating the expansion at^ 1 ^ n = 2 and |q| =2 
because r* is very large compared to R and e* is small for both the Sun 
and the Moon (other bodies do not matter here), (57) becomes (see Kaula, 
1964) 


V*(x) = 


•ii— V (JL_\ — ■ ■ ■ 

a * m=0 a * (1+6 m0 ) (2+m)! 


2 2 

P nm (sin ♦> I / 2 £I*> l O p) 

P=0 r q=-2 ■ 


cos [ ( 2-2 p+q ) ( u>*+M* )-q W*+m( ft*- 9- A) ( l-(-l) 2_m ) ] . 

(59) 

Traditionally, tidal theory has been formulated in ecliptical rather 
than equatorial coordinates, because in this system the inclinations of 
the disturbing bodies are very stable, which simplifies the mathematical 
treatment of their motions. The transformation of (59) to the ecliptical 
system is explained in (Marsh et al. , 1983). If s^g^y and s i(Moon) 
stand for the various ecliptical Keplerian elements of the Sun and the 
Moon, respectively, then I(g un ) and are both zero, by definition, 

(^To study some very minor tides, degrees 3 and 4 are included. 


and the origin of ft is the direction of the vernal equinox, as in the 
equatorial system. Instead of Keplerian elements, it is common practice, 
in tidal work, to use the following linear combinations of them: 


L * e+X-s-ir/2 


s = (OH-Mfft) 


(Moon) * 


h » (UH-M) 


(Sun) 


the mean local lunar time (expressed as an 
angle); 

the lunar mean longitude; 
the solar mean longitude; 


^ ~ ^(Moon) 


the mean longitude of the lunar ascending 
node with a minus sign; 


P = W (Moon) 


the mean longitude of the lunar perigee; 



0) 

(Sun) 


the mean longitude of the solar perigee. 


Here the word "longitude" refers to the "dog-legged" astronomical 
longitudes (ft+aH-M), etc., which are very close to the ordinary longitudes 
in the ecliptical system because of the small inclinations of the 
orbits. Like the mean arguments w, M and ft of artificial satellites, the 
mean variables s, h, N, p 0 and p^ change almost linearly over long 
intervals, in this case of the order of a century (thus the expression, 
"secular variations") so they are represented quite accurately by cublcs, 


like s * s Q + S|t + s 2 t ^ ^ Sjt^, etc. (see The "Explanations" to the 
American Ephemeris and Nautical Almanac) where the "rate" Sj is orders of 
magnitude larger than S2 and S3. The departures from linear change 
represented by S2 and S3 are mostly due to the gravitational forces 
exerted on the Earth and the Moon by the other planets. All the variables 

defined above, including the Earth's Greenwich angle 0, have mean speeds 

or frequencies that remain virtually constant over decades. 

The Keplerian elements a*, e* and I* are nearly constant in the 
ecliptical system, the Sun-Earth-Moon system being a very stable one, 
while to* , M* and fi* change at approximately constant rates, as in the 
orbit of a satellite. This is particularly true of the Sun; in the case 

of the Moon, Whose mass is much smaller, the pull of solar gravitation 

alters the orbit to the extent that, in the same way as the temporal 

variations (in inertial space) of the non-zonal part of the Earth's 

spinning field requires the introduction of the angular velocity 6 

(through -m0) in the theory of the motion of a spacecraft, the mean 
angular velocities of the "orbit'' of the "moving" Sun must be included in 
the formulation of the motion of the Moon (Brown, 1905, Brouwer and 
Clemence , ib. 1966, Ch. 12) . So the expansion cf V^ oon ^ must contain 
both lunar and solar frequencies, at least when the perturbations of the 
lunar elements are formulated using a first order analytical theory; only 
solar frequencies appear in V T ^ Sun ^. For this reason some tides are known 
as "lunisolar" , others as "lunar" , and yet others as "solar" , on account 
of the origins of their frequencies. Instead of analytical 
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approximations, one could use direct estimates of the osculating Keplerian 
elements; these elements are available (in ecliptical coordinates) from 
precise astronomical ephemerides based on long series of observations (see 
Seidelmann, 1982). This last approach is adopted sometimes for the very 
accurate calculation of satellite orbits (Williamson and Christodoulidis , 
private communications). 

Regardless of the body of origin, r*, and X* in (58) can be 
expanded in trigonometric functions of the ecliptical variables 
L,s,h,p,N,p' (Brown, ib, 1905). Considering only the second harmonic of 
V*, as in (59), and replacing and X* with their expansions, the 
result has the general form at the Earth^s surface (|xJ*R): 


yx> 


= £ P (sin 4>) I 

L nm ' T.f U 


m=0 8j=m 8 2 B 3 3 4 3 5 f 6 


WiVs?*- 


cos IP'L + $>*3 + 8 3 h + $'p + $^N + $'p ' + (m-n) -j] , 


8'^ are related to the 6^ as follows: 


8 36 8 ^ 
p i i 


( 60 ) 


3 i = K + 5 • 


(1=2,3, 4, 5, 6) 
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where the integers 


The 3^ are always positive or zero, except in some unimportant 

cases where the amplitude V is very small; moreover, 0. * m. The 

B 1 * * • % 1 
number D whose digits are the 3^: 


D * WMh 


is the Pood son number , which identifies uniquely each line in the tidal 
spectrum, or tidal constituent , besides providing complete information 
(once the various mean rates are known) on the value of its frequency 
and on the phase angle x T : 


l ’3 1 ,..3 6 ® + H s + 0 3** + H? + H ^ + ^ 


(61a) 


= ( 1 ) 


( Q + a;s ‘K erh +■ e;p + g'N + k?' 

v 3. * • . g, 1 o 2 o 3 o 4 o 5 o 6 o 

ID 


(61b) 


= X„ 


Because of the choice of L as a variable, the speed of the Moon orbit 

seems to appear in all tidal frequencies, but it cancels out with for 

the "solar tides V. The number g. g 0 g 0 g, g c g> and the amplitude V 0 0 

1 ^ J 4 j O Pj • • • 

Were computed and tabulated by Doodson (1921) for virtually all the 
important tidal constituents. The first digit 3j, which is also "m” in 
(58), ranges between 0 and 2 according to (59). Because it multiplies 
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6, which is part; of L, and 9 is much larger than the other rates 
involved, 6^ represents the separation of the tidal spectrum into three 
wide-apart main bands or species : long period (-.8 * m = 0), diurnal 

( 1) and semi-diurnal (g^=» 2). In turn, each species is divided 
into finer bands or groups , separated by intervals of ( or, 
approximately, by multiples of one cycle per lunar month. Within each 
band, the individual lines are separated by multiples of one cycle every 
8.85 years (0£p), 18.61 years ($£N)> and 25800 years B 6 'p'), (1) 

respectively, Usually, it is sufficient to add together all lines 
separated by less than one cycle per year from each other, and to consider 
these sums as the actual tidal constituents; if great precision over very 
long periods is needed, the hyperfine structure of the spectrum must be 
considered in full detail. The dominant constituents are far larger than 
the rest, and correspond to the principal tides M2 , S2 , K1 , 01, etc., in 
the classical notation introduced by Darwin in the 1880's. Some of them, 
together with their Doodson numbers, are shown in Table 3. The 
equilibrium amplitude g~* o P n ^(sin<f>) (where g « p/R^) is the 

size of the oscillation in the equilibrium figure of a perfectly fluid 
Earth that the tidal constituent would produce on its own (sometimes this 
is known as the "geoid tide"), and varies with latitude. Because of the 
continental boundaries, the submarine topography and the rotation of the 
Earth, ocean tides with periods of less than several years are not "in 
equilibrium", and have rather complex spatial structures. 

^^p' is, in fact, the rate of precession of the Earth's axis, treated 
here as a movement of the Sun's "perigee" because of the way in which the 
ecliptical coordinates are defined (i?(sun) a ^^ r 
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TABLE 3 


Principal Constituents of the Tidal Potential^ 3 ^ 
Caused by Both the Sun and the Moon 


Peak Amplitude 
of the 
Equilibrium 


Darwin's 

Name 

Doodson 

Number 

Period 

(hours) 

Tide 

(meters) 

m 2 

255.555 

12.42 

6.2423 

s 2 

273.555 

12.00 

0.1128 

n 2 

245.655 

12.66 

0.0464 

k 2 

275.555 

11.97 

0.0307 

K 1 

165.555 

23.93 

0.1416 

°1 

145.555 

25.82 

0.1005 

Pi 

163.555 

24.07 

0.0468 

Ql 

135.655 

24.86 

0.0193 

M o 

055.555* 

CD 

0.1356 

S o 

055.555* 

00 

0.0629 

Mm 

065.455 

661.31 

0.0222 


fne periods are those of the main constituents of the semi- 
diurnal, diurnal and long period bands. ’’Equilibrium tides" of 
semi-diurnal period peak at the equator, diurnal ones, at mean 
latitudes, and long period ones, at the poles, 

"Frozen" lunar and solar tides (respectively). 
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In what follows, expression (60) will be written in the following 


simplified way: 


V^(x) - l V T P 2m ( sin4>) cos(o) T t + X T +inX + (tt-n) ~)] , 


(62a) 


where indicates the sum of all significant tidal constituents, and"*" 
T 

has been dropped from the right hand side, as the solar or lunar origin 
of a particular constituent is of no further interest here. Now, 
according to (61a) and to the definition of L, is a function of ur; 

m ■ Int [&/|w | + 1/4] , "Int [ .]" meaning "the integer part of". In the 

T 

A 

new expression of Vj(x) , 


T e i S 2 S 3 e 4 e 5 e 6 


(62b) 


(Pj * m) 


is the amplitude , 


°T * 3j(6-s) + + + £|p + N + fi'p' 


(62c) 
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is the (angular) frequency and 


M 0 -s -*/2) + Kb + Kh n + K? n + O + Kp' 
loo lo 3 o 4 o 5o 00 


(6 2d) 


is the a stronomical phase . It Is a useful and universal convention 

that X T is the phase at t' * 0.00 hours , Greenwich mean time, of the day 

in which the instant "t" falls (so "t", in fact, is t-t**). But in the 

o 

discussion that follows, in order to make the treatment of tides 
consistent with that of the radial orbit error, X t is defined as 


X T * xJt .) 

(as in (62d)) 


- u) t 

TO 



(62d ' ) 


where ’’tg" is t Q for the general case corresponding to (47), but (as 
explained in the comments that follow (52)) for "frozen", repeat orbits t Q 
is such that M(t Q ) ■ 0 at the start of the "repeat" part of the mission. 
This choice is preferable for writing formulas: for actual calculations, 
the conventional one is more accurate and freer of numerical round-off 
problems. 
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The change in the geographical coordinates of the point where the 
spin axis intersects the Earth's surface, or polar motion, reflects mostly 
the free wobble of the Earth (Chandler Wobble). Polar motion produces a 
fluctuation in the centrifugal force experienced by objects fixed to the 
planet, so it changes the acceleration of gravity acting on them, and 
causes deformations and movements in the body, oceans and atmosphere 
similar to those due to the tidal potential (but much smaller). This 

response is known as the pole tide. Since its period is completely 

unrelated to the lunar and solar frequencies (about 1.18 years), it can be 
detected in sensitive gravimetric tidal records as a separate and very 
weak oscillation, particularly at high latitudes. Precession and nutation 
also affect gravity, but they are themselves parts of the overall tidal 
response to the pull of the Sun and the Moon, and their effects on the 
solid Earth, oceans and atmosphere are incorporated into the diurnal tides 
(see Wahr, 1982). 

The tidal potential V*, according to (57), is rotationally 
symmetrical about the line from the Earth to the disturbing body. Its 

equipotential surfaces have an interesting geometry; the surface V* » 0 
is a double cone with its vertex at the geocenter (opening both towards 
and away from the body). This cone, which gets broader as the body get3 
closer, is the asymptotic surface of a family of revolution hyperboloids 
of two sheets, which are the other equipotentials. The tidal acceleration 
VV* has the same rotational symmetry and is also equal in magnitude and 
opposite in sign at points equidistant from the geocenter about the 
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Earth-body line. This description holds provided that terras in (57) of 
degree n > 2 can be neglected and one does not get so close to the 

disturbing body that it no longer can be treated as point-like. In the 
case of a satellite free-falling around the Earth, the field of "residual 
gravity" measured in a system fixed to the spacecraft has also this 
configuration (with the Earth as the disturbing body). Detailed and 
charming discussions of this geometry can be found in two papers (1977 and 
1982) by the late A. Marussi. At the Earth's surface, the combination of 
the geopotential V with V* creates, in the nearly spherical equipotential 
surfaces of V, a slight deformation shaped as a symmetrical ellipsoidal 
bulge pointing towards the disturbing body; this bulge Is the total 
"equilibrium tide". Since the Sun and the Moon lie at or near the 

ecliptic, this bulge is tilted with respect to the Earth s spin axis, so 
the rotation of the latter carries the points of its own surface through 
the bulge. At mid- and low-latitudes the points cut through both "lobes" 
of the bulge, experiencing a twice-daily variation in V*. Closer to the 
poles, only one lobe is cut (most of the time) so the diurnal species 

prevails far from the equator, and the twice-dlurnal closer to it. As the 
mean value of V* is not zero along parallels, there is a component that 
changes only as the body moves slowly along its orbit: this is the long 
period species. The ellipticity and precession of the orbit modulates 

both the bulge and the species, splitting these into bands and those, in 
turn, into lines. 
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It is Interesting to notice that, according to (57), the tidal 

potential increases as r* decreases. Imagine that the system of 

coordinates is chosen at the center of mass of a body which, like the 
Moon, orbits a more massive one, like the Earth, and that the two bodies 

get closer, until the attraction of the main one exceeds the self- 

gravitational force chat pulls the other together, and breaks it up. The 
distance at which this catastrophe occurs is known as the "Roche limit". 
Maybe this is the process that once created the rings of Jupiter, Saturn 
and Uranus out of the disintegration of icy moons, and now keeps those 
rings from coalescing again into larger bodies. In the case of Io, the 
innermost of the large moons of Jupiter, the tidal forces are sufficiently 
vigorous to heat the interior by continuous friction, probably keeping it 
partially molten and causing the spectacular vulcanism for which this 
world is now famous. Terrestrial tides, though much less energetic, play 
a significant role in the deceleration of the Earth-Moon system. 

5.2 The Tidal Response 

The attractions of the Sun and the Moon torque the slightly 
elliptical spinning Earth, causing the precession of its figure axis, one 
cycle of which takes about 25800 years. The relative positions of the 
Earth and the disturbing bodies vary as they move, so the torque is 
modulated at the orbital frequencies both in intensity and in direction. 
The result is a "nodding" motion of the Earth's axis, or nutation (see 
Leick, 1978). Precession and nutation are movements of the Earth as a 
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whole. Because the planet is not perfectly rigid, the changes in the 
attractions of the disturbing bodies also produce small movements, or 
tides , which vary from place to place, as well as slight variations in the 
spin rate 5 (not considered here any further). These tides are the 
subject of the rest of this section; the reader can find a very clear and 
up-to-date discussion of the (somewhat arbitrarily separated) components 
of the tidal response of the Earth: precession, nutation and tides, in 
(Wahr, ib. 1982). 

5 . 3 L inearity and Time-^Invariance 

The joint attractions of the Sun and the Moon cause displacements of 
the matter in the interior, oceans and atmosphere of the Earth that are 
very small compared to the size of the planet; that joint pull itself is 
quite small when compared to the gravitational force the whole planet 
exerts on each of its parts. It is, therefore, not surprising that those 
displacements, their interactions, and the resulting fluctuations in the 
gravitational field can all be described accurately by first order, or 
linear, approximations to the response of a nonlinear system (the Earth) 
to the disturbing forces. Over an interval of a few centuries, the large 
scale mechanical characteristics of the Earth are likely to stay virtually 
constant, so it seems reasonable to expect those linear approximations to 
be also time-invariant , perhaps after making some minor simplifying 
assumptions. If such is the case, the spectra of the various components 
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of the tidal response (ocean tides, body tides, etc.) should have thej r 
lines situated at precisely the same frequencies as their driving forces, 
which are those given by expression (62c). This is in very close 
agreement with observation (except for some small anomalies in shallow 
seas that suggest some nonlinear frequency mixing) , which means also that 
both the disturbing forces and the frequencies of the resulting tides can 
be calculated very accurately from purely astronomical data. 

The formulation of a realistic nonlinear model for the Earth, and its 
linearization with respect to some reasonable and convenient undisturbed 
state to obtain a good time-invariant approximation, are quite complex 
matters, as is the solution of the resulting linear partial differential 
equations. Nevertheless, because of their time-invariance, and of the 
nearly spherical symmetry of the Earth, the linearized equations have 
solutions with strikingly simple and helpful properties, summed up in the 
ideas of "admittance function" and "Love numbers". 

5 . 4 The Admittance Function 

A long record of observations of the "surface-bottom" ocean 
tide, at a station of latitude and longitude (f, A), can be approximated 
quite accurately by a trigonometric series: 
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- l l a s -CM) cos(t» t + X. - B 7) 
0 Tp=0 1 c 


“ 1 c (4»,M) , 

T 


(63) 


where W T is the tidal frequency and X T the astronomical argument* as per 
(62 c-d) , of the ocean tidal component ? qT » Now, according to (62a) 


V TT (M»t) = V T P 2m (sin 4») cos [m T t + X T + tn* + (m-n) 


(with m = Int f 5/ 1 a> T | +4-]) 


(64) 


is the constituent of frequency of the tidal potential at ( <J> , X ) . Over 
the Earth's surface V^, T consists only of second degree spherical 
harmonics, under the approximations made in its derivation. The 
admittance function relating to V^, T is the complex quantity 


Y(<M,« ) = 


U oX (<M)l ie T (i,X) 


I,v Tt (M)I 


(65a) 


where 


- tM a $ T ( *» A)jl/2 si ® n f v r ] 


(65h) 


3=0 
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and 


e T (<f>>^) « -tg 1 [a lT (<j/,X)/a 0T (<j) f A)] - mA •» (n-m) j 


(65c) 


are the amplitude of the constituent of of frequency and its 

phase shift relative to the corresponding constituent of V*. The 
accumulated evidence (Hendershott and Munk, ib. , 1970) suggests that 
Y(<}>, A.u^) is a smooth function of w , so it varies very little in 

modulus and phase across each of the narrow bands that muke up the tidal 
spectrum. For this reason, if ai T ( centra i) is t,1e frequency of a line 
within one of these bands situated near the center, it is valid to assume 
that 


S<M) - S( C e„tral> ( + ,W 


(66a) 


and 


C (<J>,A) 

OT 

V (<j>,A) 
Tt 


S <<f>,A) 

oc(central) 

V ( <l>, A) 

Tt( central) 


(66b) 


so 


a (if), A) a (<j>,A) 

0t Ox(central) 

a (<f>, A) a (I, A) 

It 1 t( central) 

for 0 < 3 <1 and all w in that band. 

(63) follows that 


(67) 


From the last relationship and 


5 ■ i a (M> l (V /v , n ) cos(O) + X T - 3 f) » 

o(band) 3=0 3x(central) x(band) T(central) 2 


where "l" indicates that the sum is over the constituents in the band 
t(band) 

only. To speed up calculations at the price of a slight loss of accuracy, 
one can make use of the relationships 


cos(<*> T t + x T - 3 j) = cos x T cos(0) T t « 6 j) - sin X T sin(w T t - 8 j) 


and, to first order in (w - to ) t , 

T( central t 


^kV-s!) - :^k^ (ce „ trant -^H(S(central)-^) c Cs!<“,(ce». ra l) t - 3 5)- 


So, instead of (68), 


’o(band) 


C<M>t) - l a (4,X) {A (t) cos(w i 

3=0 3t(central) x(band) x(central) 


+ B (t) sin (io 
x(band) 


x(central) 


- 3 - 1 )} , 
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where 


A (t) 
x(band) 


! 

t(band) 


( VV*«»twl> tfl *t<«.ntMlf V* sin *t + cos V 


<70a) 


and 


B (t) 
t(band) 


l (v/v 


.)[<« 


t(band) 


T r(central) r(central) t 


- W_)t cos X - Sin xJ 


(70b) 


provided that l(» T(c „ tral , - “ t )t| « K“ l(ce „ tral) H X,|. In actual 

calculations this condition can be met by updating the phases X T In the 

above expressions sufficiently often (for a given band), and using these 
new values with "t" counting from the time of the last resetting, until 
the next update. This saves the effort of having to compute one sine and 
one cosine value at every instant where (6b) has to be evaluated, for 
every frequency in the band,. For solar tides, the V must be corrected to 
account for the "radiational tide potential” (Cartwright and Tayler 
(1971)), Radiational tides are caused mostly by the periodical heating of 
the atmosphere by the Sun, which modifies the air pressure at the sea 

surface, forcing the ocean at the same frequencies as the solar tides. 

This thermally driven variation of the atmosphere accounts for most of the 
atmospheric "tides", with the gravitation of the Sun and the Moon playing 
a lesser role. 
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5.5 Love Numbers 


5 can be regarded as a function defined for all <j> and X but which 
o 

is identically zero over land. Such a function must have an infinite 

spherical harmonic expansion because, although very smooth over the 

oceans, it is discontinuous across the coastal boundaries. Furthermore, 

£ is at each point a sum of trigonometric functions of time, as in (63), 
o 

so its complete expansion in space and time is 


UM, t) - l | { l I l C nma g T Y nmct <M)cos(to T t + X - " $ T> 1 


t n^O m=0 et=0 3=0 


* \ \ t onT ( *- X > t) > 

t n=u 


where the 


5 o are the spherical harmonic Coefficients of a„ ( <j>, X) 
nmotpt pT 


(3=0,1) in (63).^^ Putting (70) and (71) together: 


^o(band) ( X * C> “ J 0 J Q Jq 3=0 ? ™a3x( central ) Y nm e ( <M) 


{A ( t) cos(a) n a . t - 3 j) + B (t) sin(u > .t - 3 ■?)}. 

t(band) t( central) 2 T(band) T(central) 2 


(^Tide-lands not withstanding, in a global scale it is reasonable to 
regard the coasts, as well as the ocean floor, as impermeable and well 
defined boundaries. Under this assumption, the conservation of mass 
requires that £ n g=0 in (71) et seq. , so foOQ/Jr for all ft and all T 

(if 0> r /0). However, the same notation shall be used later for the 

spherical harmonic expansions of tidal models and their errors; some of 
them do not satisfy the conservation constraint, so the zero harmonic has 
been- kept in (71) et seq. 
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The altimeter measures the total tide, not individual components, 

and the concept of admittance makes it easier to take into account the 

’’sidebands" when estimating the "main tides" from altimetry (such as M2, 

S2, Kl, etc.). In what follows, 5 ( 4>, X, t ) (the expression between 

oni 

curly brackets in (71)) shall be called, for short, the "n spherical 
harmonic" of (at frequency ur). 


The tidal displacement of water changes the forces acting on the body 

of the Earth, which yields because this body is not rigid. These related 

movements of water and solid matter modify the gravitational potential. 

Since the responses are nearly linear and the Earth is close to having 

spherical symmetry, the spherical harmonic expansion of is related to 

those of the vertical displacement 6 ' of the bottom and of the tidal- 

(? o ) 

induced change in potential V as follows: 


6^(4>,X,t) = 1 l a h* 5 6nT (i,l,t) 
t n=0 


(73) 


and 


(?) 00 n+l 

v ° ( 4 >,x,t) - g I l ) 5 onT (M,t) (f) , 

t n=0 -v 


(74) 
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Here 


3p 

w 

a * , 

n p g (2n+l) 

2 

g »■ p/R is the mean acceleration of gravity on the Earth's surface, p 

W 

the mean density of water and p the mean density of the Earth. h“ 

e nt 

and k'„ are known as load Love numbers (Munk and McDonald, 1960; 

. . nx . . . . ■ 

(V) 

Farrell, 1972) because 6 and V ° are the result of the loading of 

the body of the planet by the ocean tide on its surface. Their dependence 

on to' as well as on n is stressed here to bring the formulation in line 
with that of Wahr (ib,, 1979) for ordinary Love numbers, which probably is 
the most realistic one at present; the extension to load Love numbers can 

be found in (Sasao and Wahr, 1981). In the past this dependence has been 

largely ignored, because the theory was restricted mostly to the 
spherical, non-rotating case, where the main free oscillatory modes have a 
frequency considerably higher than any tidal component. Wahr has included 
the effect of both rotation and oblateness in his computations. The 
latter requires the use of functions slightly different from 

(<M) 

the Y ; however, the difference in the results is less than 1% and 

nmo 

can be ignored here. 

The tidal potential has an expansion on the Earth's surface given by 
(62a). V* causes a deformation of the body of the Earth, or body tide, 
resulting in a vertical displacement 6 and a change in potential V' 
which are, approximately. 


(here v Tt2 < X, t ) is the same as V Tt (<}>,M) in (64)), 
and 

vSIrM) - I h 2T V Tx2 (4>,>v,t) (76) 

. t 

The contribution of C Q to the geocentric tide of expression (1) is 

5o - 4 - «o + 5 ' 

m l \ (1 * “n h „V W*’ 4 -'* • (77 > 

t n=0 

The h 2t and k 2x are known simply as Love numbers , after the man who 
first wrote about them (Love, 1909). In the interval 2 < n < 25, which 
includes most spatial frequencies of interest, the approximate sizes of 
these various numbers change monotonically from 

h 2x « 0.612, k 2t * 0.302, h' T = -1.007, k' x a -0.310, « 2 a 0.11, 

to 

h 25t “ 0-047, k 25t a 0.002, h^ T a -2.194, =-0.046, = 0.011 

(from Table 1 in Hendershott, 1972). The largest departures from these 
values occur within the diurnal band. 
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The values of k^, h^ and k£* W can be calculated by using a 
simplified linear model for the departures of the particles that made up 
the body of the Earth from their undisturbed state, selecting a plausible 
set of physical parameters for the matter of the planetary interior, and 
integrating numerically a set of ordinary differential equations related 
to the linear model. The integration variable is depth, in the interval 
from the geocenter to the surface, and the Love numbers are the integrated 
values at the surface. Usual approximations include a spherical or 
elliptical body where density, elastic parameters, etc. vary only with 
depth and viscosity is zero, so there is no energy dissipated by the tidal 
motions. The mechanical parameters are discontinuous at the boundaries 
separating the core (divided in solid inner core and liquid outer core), 
the mantle and the crust, sometimes with an additional change between 
upper and lower mantle. In consequence, the published values of the Love 
numbers vary according to the linear Earth models used for their 
computation. Wahr (ib., .1979) has assumed an elliptical, vertically 
stratified, dissipationless Earth with no oceans. A slightly more 
accurate formulation, including dissipation, would require the 
introduction of small additional phase shifts on top of the X T in 

^ The re is a component Vp of the total tidal potential, in the long 
period species (ra=0) , which is actually constant in time (and virtually a 
second zonal in space, see expression (62)). It produces an equally 
constant change in the Earth's equatorial bulge and, thus, in the value of 
Cgoo compared to what it would be if our planet were alone in space. The 
Love numbers for the resulting stationary body tide are probably those of 
a "liquid Earth", as our planet is supposed to have been "flowing" towards 
its actual figure of hydrostatic equilibrium over the eons that this 
"frozen tide" has been acting on it. There is also a slight discrepancy 
between the "bulge" of the mean ocean surface and that of the geoid , equal 
to the "frozen" component of the ocean tide Vp/g (S 0 and M 0 in Table 3). 
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the arguments of the tidal functions. This appears unnecessary here. The 
main variation of the Love numbers with o)^ occurs when the frequency is 
close to one cycle per day (diurnal tides) and is caused by a resonance in 
the fluid core (Wahr, ib., 1979 and Wahr and Sasao, 1981). For semi- 

diurnal tides the Love numbers are virtually independent of ur. 
Horizontal changes in density, etc., further complicate the picture, but 
they do not have to be considered here; they would only matter in highly 
accurate or local studies of tides. 

5.6 Representing The Ocean Tides 

The mathematical model for the ocean tides is derived from the 

nonlinear Navier-Stoke's equations of hydrodynamics, linearized by 
ignoring small quadratic terms in the unknowns and using first order 
formulations for the forces related to turbulence, bottom friction, 

crustal bending (6") and self gravitation (V^^) . In Cartwright 
(lb., 1977), Hendershott (1972), Schwiderski (ib. (1980b)) and Parke and 
Hendershott (1980) one can find information on the historical development 
of the linear theory and of the methods for integrating the tidal 
equations.' Alternatively, the use of nonlinear terms to represent bottom 
friction yields a system of equations that is overall nonlinear and can 
only be solved by numerical integration in space and in time (Estes, 
1980). On the other hand, a purely linear system requires numerical 

integration in space alone, the time variable being eliminated by working 
with the Fourier transform of (as in Pekerls and Accad , 1969). The 
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simplest linear equations are those of Laplace. Usual assumptions made in 
obtaining tidal equations are that the Earth is spherical, the depth z of 
the ocean is nowhere significant compared to the planetary radius and to 

the horizontal size of the tidal waves, that the small meridional 

component of the Coriolis force can be ignored, and that the velocity of 
the horizontal tidal current is independent of depth. In fact, the 

velocity does change from the surface down, but because the density of 
ocean water is well stratified there are two main types of waves: those 

described by the tidal equations, which move the surface up and down to 
produce C Q » creating horizontal currents that are virtually independent of 
depth, and those that do not move the surface appreciably (also known as 
internal waves) but cause currents that vary with depth. Clearly, only 
waves of the first type can be mapped with an altimeter or with 

conventional tidal gauges. 

Figures 5 and 6 show the charts of the M2 constituent obtained by 
Schwiderski (1979). The lines in Figures 5 and 6 present the amplitude 
5 ot O,A) in cm and the phase in degrees. The lines of equal phase (or 
cotidal lines) in Figure 6 show, successively, where the tidal waves crest 
at about half an hour intervals. Looking at these lines in their proper 
Sequence makes an ‘’animated picture" of the waves as they turn about their 
fixed nodal points or amphidromes (where the amplitude is always zero). 
Notice how smooth the Waves are. 
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Goad (1980) has computed the power spectrum of Schwiderski's M2 . tide 
up to degree n-180, showing a fast rise up to degree n-8 and then a fall 
of more than an order of magnitude at n-20; from there on the spectrum 
trails off to zero slowly, reflecting the discontinuities of along the 
coasts. For the other semi-diurnal components the general shape of the 
maps, including the approximate location of the amphidromes , is quite 
similar as for M2. The same is true for the diurnal components, although 
their patterns are quite different concerning the locations of the 
amphidromes, etc., from those of the semi-diurnals. 

While the waves of individual constituents are virtually periodical, 
the total tide is not, because the orbital frequencies of the Moon are not 
harmonics of those of the Sun. .therefore never repeats itself exactly. 

A fluid motion can be described in terms of a velocity potential <j> 
and a stream function <|> (see Lamb (1932)), and both and can be 

expanded into what can be a fast converging series of velocity potential 
base functions <f>^ (Proudmatf, 1916; Rao and Schwab, 1976): 

{.«?}($, l,t) “II I { J£ T } MM^osU + x, 

? 0 T i=0 0-0 -igT 

where the and A ^igt a,re rea1- nura,5ers • 


- 6 f) , (78) 
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The are specific to a given ocean basin* as they depend on the 

bathymetry and the coastal outline, and form a complete set of orthogonal 
functions over that basin, in the sense that 

//(basin) \ d0 " D i • (79 > 

where the positive constants are zero if i f j, 
and 

//(basin) *i 0 < 80 > 

for at least one value of i, unless ? is identically zero. Integration 

0 

here is confined to the oceanic regions, and do is the spherical element 
of area 

do « cos<t> d$ dX . (81) 

Notice that the are independent of u> T ; in fact, they can be 

used to describe a wide variety of vertical motions of the sea surface 
within a given ocean, besides tides. These functions do not have 
analytical closed expressions except when the shape of the ocean is 
unrealistically simple. To compute 4^(4, X) one must solve numerically an 
eigenvalue-eigenfunction problem involving a differential operator related 
to the unforced linearized tidal equations, with the boundary condition 
that, along the coast, 
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0 


(82) 



Here, "z" is the depth of the ocean, and "3/9n" represents partial 
differentiation in the direction normal to the coast. The calculations 
can be done only with a digital computer, but, once the have been 

obtained, this disadvantage can be offset by the speed with which the 
expansion (78) of t converges and by several benefits associated with 
orthogonality. Rao and Schwab (ib. , 1976) and Sanchez et al, (1984) have 
Studied the use of these interesting functions for parameterizing tides 
and other oscillatory motions of large bodies of open water without having 
to make assumptions on the amount of internal dissipation and bottom 
friction, usually poorly known. The 4^ can represent any long-wave 

vertical displacements of water, even when the energy dissipation is 
unknown, in much the same way as spherical harmonics can be used to 
represent any gravitational potential, even when there is no knowledge of 
the mass distribution that generates the field. 

If the bottom of the ocean is replaced with a horizontal surface 
lying at a depth of "no motion", the 4^ calculated with this new boundary 
can be used to parameterize the stationary sea surface topography w^ 
corresponding to the global circulation. 
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5.7 The Tidal Part of the Orbit Error 


The attractions of the Sun and the Moon, together with the tidal 

changes they induce in the gravitational field of the Earth (expressions 

(74) and (76)), affect the orbit of a satellite and must be considered 

when computing its position and velocity. The direct influence of the 

celestial bodies can be calculated with such accuracy that the results 

can be regarded as exact. The solid Earth and ocean tides are less well 

known, so there may be significant errors in the calculated values 

of V( 5 ) and . The body tide 6 has a limited frequency range, mostly 
o 

below two cycles per revolution, because only the second harmonic V* has 

an important effect (expression (76)). As the existing models of 6 are 

quite accurate, only the long period effects of their small errors AS are 

likely to matter, adding to the polynomial terms and to the "quadratic" 

oscillations in Ar (expressions (47) and (52)). Therefore, shorter 

period errors related to tides are probably due mostly to A? , the 

o 

uncertainty in the ocean tidal charts. The perturbations of spacecraft 
orbits have been used to learn about the broader features of tides, and to 
try to estimate directly Love numbers and other important geophysical 
parameters (Lambeck et al. (1974), Felsentreger et al. (1979), Marsh et 
al. ib. (1983)). 


The error A? can be represented by an expansion like (71) for C , 
o o 


With the coefficients 


t, „ substituted by the errors 
nmaPx 


Ac . in 

nma0x 


these coefficients. Replacing with its equivalent in Keplerian 

nmot 
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elements (expressions (29-30)) and using the relationships for the 
/roducts of sines and cosines. 


~ 2 l ^ fi ^ n AC nmagT ^ F nmp ^ 

t n=0 m~0 a=0 3=0 p=-n 


cos [( (n-2p) (aH-M) + m 6 ' ± o) T )t + 4> nm ap0o * (X T "3 7 )] * 

(83) 


where apOo 33 C n “2p) ( a>(t:^) + M(t^) ) + mO'Xt^) ; the osculating w, M and 
Q' have been approximated by &t + u)(t ), etc,; and X T i» as in (62d^). 
Because of the nearly linear relationships between tidal phenomena, 
rewriting expression (74) with the harmonics of A? Q instead of those 


of C Q gives the effect AV^ A ^°^ of AC Q along the orbit. So, 
replacing (l ) n+1 Y nma (<j>,X) in (74) with its equivalent in Keplerian 
elements according to (31), gives (Lambeck et al. , ib. , 1974) : 


( A ? ) „ 

I « 0 o*k' T ) (|) 

2R nmotSt w 


n +1 


AC 


nmaBT ~ nmp 
p=-n 


l F (I) l G 


(e) 

npq 


cos [ ( ( n- 2 p+q ) ( uH-M) - qo) + m9' ± oj t ) t+<f> nraapqo ± (x x _ 3 %)] . 

(84) 
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Neglecting the contribution of the error in <$, the sum of the shorter 
period perturbations that form the "periodical" part Ar^ of the tidal 
radial error is given by 



(t) - 


_ 1 _ 

2R 


l 

nmpqot3x 


°n ()+k nT )A; nmaf5T ^ r nmp(q+{±l})T 


cos[((n-2p+q)(uH-M) - qu) + m9" ± U))t + 4’ nmap ( q .p{. fc i }) a ± (X x "3 7 )]} * 

(85) 

where "(q+{±l}") stands for "q" , "(q+1)" and (q-1)", respectively, and 
"(T)" for "tidal". The frequencies present in Ar^ differ from those 
in Ar v (expressions (47) or (52)) by ± u> T , so the actual values of 

the coefficients r nra p(q±[ 1 } ) oX are a ^ so different (see (16), (38) and 

(43)), thus the subscripts "nmpqx". Terms where the are substracted 

are prograde ; the others are retrograde (the change in the argument has 
the same or opposite sense to that caused by the rotation of the Earth). 

(NT) 

As in expression (47) for Ar v ' , the mean elements and their rates 
generally change, very slowly, as the orbit decays, though they may be 
considered as constants during each weekly arc; the r n mp(q+{±l} ) change 
accordingly in (85) . 
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5,8 The Case of the "Frozen”, Repeat Orbit 


When o> « 0 and exactly revolutions take place In precisely Np 
nodal days, the ground track repeats itself with a frequency 
w = 2'ir/T j where T rc “ x (length of the nodal day). Then the 

frequencies in the arguments of the cosines in (85) are, with the possible 

/ T ) 

exception of w , multiples of u) r so Ar becomes 
T rc 




2 R nmpqaST(nonres) + * tRT ^ Wq+I±° }> T 


cos[(((n-2p+(q+{±l}))N R -m N D )u rc ±0) T )t+<|> ninap(q+ | ± 5 j )o ± x^gl)], 

( 86 ) 

When dealing with this type of orbit, the time origin t Q for the 

^nmap(q±{l})o and x t is as in expression (52) and in (62'c) , 

so M (t Q ) =0 early in the repeat part of the mission. 


The error At in the geocentric tide is present in the residual 
altimetric saa surface heights according to (7). Assuming that the error 
in the body tide is small enough compared to A? q to be neglected, 
expressions (29) , (30) and (77) imply that 


4 ' G<t) “ 2 Jist pl-w (I) 


cos[((n-2p)(oH-M) +m8" ± w T )t + 


± Cx-ef)]. 


In the case of the "frozen", repeat orbit this formula becomes 


" * „ie t pi/™, <« 


(87) 


cos [«(n-2p)N R - m tf p j ^ ^ t + 4> nmap0o * (X T “3 f ) 1 • 


( 88 ) 


Comparing (86) and (88) one sees that the periodical terms in both 

~(T) 

Ar and A have the same frequencies when these errors are treated 

as time series. Consequently, Af(T) and A?^ are "lumped together" in 

the residual sea heights and can be separated only by using the Love 

numbers kf- , so these must be known reasonably well. The same comment 
nx 

can be made regarding the h*^ needed, in turn, to obtain AC Q from A?q, 
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5.9 The lasing of Ocean Ties With the Mean Sea Surface and With 
Each Ocher 


When observed with an altimeter from a "frozen", repeat orbit, a 
tidal constituent of frequency <*> T would appear, at the same points along 
the repeating ground track, to rise and fall between consecutive 

overflights at the slower frequency 


U) 


t( st robed) 


* 0 ) - at Int [w /u> + 1/2] . 

t rc t rc 


(89) 


(Int [x] = integer part of x). This phenomenon is analogous to the 

apparent slowing down of a rotating wheel, or the plate of a record 
player, when illuminated at regular Intervals by short flashes of light. 

The "strobed" frequency W x( strobed) wou -*- £ * zel|:o if were an exact 

multiple of the repeat frequency o> .. The corresponding tidal 
constituent would appear "frozen" in time and inseparable from static 

features of the sea surface such as the geoid undulations or the 
stationary sea surface topography. This may happen with some tides 
depending on the inclination, eccentricity and height of the orbit 
(a) = (0-ft)/ N n , and ft depends on ~a, "e and I according to (10)). 

Moreover, two different tides of frequencies a> T ^ and may appear 

"lumped" together if their "strobed" frequencies are the same. For 
small "e and altitudes of about 1000 km this is possible when the 
inclination is somewhere between 70° and 110°, which is one reason why the 
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projected orbit of TOPEX has been chosen with an I ~ 64°. In the case of 
SEASAT, where I ~ 108°, some tides (like 01 and N2) were aliased with 
each other and some (like PI) with the mean sea surface, as mentioned by 
Mazzega (1984). Details of this problem are discussed at some length In 
two reports: one by NASA's TOPEX Science Working Group (1981, Appendix 
B.2), and the other by a similar group organized by ONES in France to 
discuss POSEIDON (lb., 1983). 


5.10 The Complete Expression of the Radial Orbit Error 


The total error Ar is the sum of the tidal and non-tidal components 


Ar^^ and Ar u,A-> (expressions (47) and (85)) 5 


Ar(t) »' Ar 




aXq(nont eS ) 4Cnm0 ‘ 


cos [ ( (n-2p+(q {*!}))( aH-M) - (q+{±l}) *a»fm*0)t + <l> n , n ap(q+ {*1 }> 0 ^ 


+ 2R nrna0qT(nonres) “" <1+k ^ )i5 nm>»6T %»p<q+{«i> T 


cos [ ( (n-2p+( q+{±l }))( wfM)-( q+ {±1 } ) «+m 9 0 1 ± to T t 


+ ^nmap(q+{±l })o ± ( V 0 f )J + V 08 
+ B^sin Mt + cos Mt + D^t sin Mt 


Mt 


+ E. cos Mt + Fi-t^sin Mt + V r., (t-c k ) , 
k k J=0 ° 


(90) 


where the terms that modulate the sin Mt and the cos Mt, and those in 
the polynomial at the end, are the agregate of all the very long period 
and resonant contributions from errors in the models of both gravitational 
and non-gravitational forces, as well as uncertainties in the initial 
state. Expression (90) describes all significant elements of Ar over a 
weekly arc. 


The non-tidal part of the error In (90) has the same structure as a 
sum of additional tidal components of frequencies 0, w, 2u>,. . . , 
(corresponding to q=0, J. ,2 , , , . , with their amplitudes becoming negligible 
above 2w), one of which is "static" (q-0). 


The corresponding expression for the especial case of the "frozen", 
repeat orbit is, from (52) and (86), 


Ar(t) 


nmctpq (nonces ) r nmp(q+{*?|) 

cosU(n-2p+(q+{±l l^)^R +mN D ^ w rc t+ ^ nma p(q+ 

a (1+k' )AC 
n v nt ' -n 

cos [ ( ( ( n-2p+( q+ {±1 } ) ) N^+ mNo) u> rc ± w T )t 


+ 2R nmaPqT(nonres) V *’'V ' 05 n»«8T r „mpCl+(±ll)T 


+ ^nmap(q+{±l })o * ^ T ~8 2 )^ 

+ A^cos Mt + sin Mt + C^t cos Mt 


+ t sin Mt + E^t cos Mt + F^t^ sin Mt 

i lf j 

+ I v k Ct~t l <) . 

J=0 JK 0 


(91) 


Here, because the orbit is stabilized by occasional maneuvers, the 
r nmp(q+ {±1 }) ant * Ml* can be regarded as constants over the whole 
mission. 


6.0 IMPLICATIONS FOR THE ANALYSIS OF ALTIMETRY 

6.1 General 

Expression (7) shows that the residual sea height Al 1 contains, 
besides the orbit error (Ar), information on the unknown parts of the 
geoid (AN) and the tides ( A^), the non-tidal variations (v t ) and the 
mean sea surface topography (w q ) of the general circulation* 

The fine details of the geoid can be determined very clearly by 

mapping the mean sea surface with altimetry, because the orbit error Ar 
and the topography w q introduce distortions at much longer wavelengths. 

Those details often outline quite well the submerged trenches, ridges and 
mountains that produce them because of incomplete isostatic compensation. 
This was recognized soon after the first use of an altimeter in SKYLAB 
(Leitao and McGoogan, 1975; McGoogan et al., 1975). Today, high 

resolution maps of the sea surface, like those produced by Marsh and 
Martin (1982), reveal aspects of the ocean crust of great interest to 

geologists and geophysicists (Watts, 1979). 

In addition to being useful over the sea for studying the 
gravitational field and oceanographic phenomena, altimetry is valuable to 
some extent over land (Brooks, 1981) and particularly over ice (Brooks et 
al, , 1978) . The main problem on land is the roughness of the surf ace 
observed; a normal sampling fate of one measurement per second (about 


every eight kilometers along track) would be Insufficient for mapping 
mountain areas, for example. Other applications, such as the study of 
wave height, surface winds, etc. (Mognard et al. , 1984), based on the 
shape and Intensity of the returning radar pulses, rather than on their 
timing, are not affected by the orbit error. 

6.2 Crossover Points and Overlapping Arcs 

The difference between two heights measured on the same spot at times 
tj and t 2 is, according to (7), 

W e 2> ■ v e i> - w 

" 4r<tj) - 4r(t 2 ) + A? G (t 2 ) - A^Ct^) + w t <tj) - w fc Ct 2 > 

+ AA(t 2 ) -■ AA(tj) . 

(92) 

The permanent part AN + w is totally absent, and only Ar, the 

o ■ ' . . 

■ ■ • •* 

tides and other time-varying features are observable to some extent. Of 
course, measurements very rarely happen to be ’taken so close to each other 
that the cancellation of AN + w q is complete. Normally it is necessary 
to interpolate Ah^tj) and Ah w (t 2 ) from measurements just preceding 
and just following passage over the point in question. There are two 
types of point: a crossover point is the intersection of an ascending and 
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a descending arc of the ground track (each named after the nearest 
equatorial crossing) „ while an overlap point lies where two or more 
parallel arcs (or colinear arcs) run on top of each other; the height 
difference in (92) shall be called a crossover difference , or an ov erlap 
difference , accordingly. In a "frozen", repeat orbit, all points are 
overlap points, but only some are also crossover points, 

Over arcs of a few thousand kilometers Ar can be represented 

closely by a constant and a linear terms, or "bias+tilt". For many years 

now the "bias+tilt" model has been used for estimating Ar from crossover 

differences, by least squares adjustment, in order to correct the data 

(Hummel and Rapp, 1977, Mather et al, , 1977). Extensions of this idea, 

involving either polynomials or Fourier series, are explained in Goad 

et al. (1980). Similar procedures for correcting overlap differences were 

also developed in the early days (Mather et al. , 1978), A system of 

observation equations for estimating biases and tilts out of differences 

of either type is rank-deficient, and to solve it is necessary to "fix" 

first some carefully selected "master arcs", or to minimize the 

discrepancies along each arc between h„ and N using some model of the 

c c 

geoid (see Rowlands, 1981). 

As explained earlier, in a "frozen", repeat orbit the radial error 
contains a significant component that is a function of latitude alone, and 
unobservable in crossover differences. Therefore, application of the 
"bias+tilt" method may leave zonal "wrinkles" in the estimated mean sea 
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surface. Fortunately, roost of the data available today were obtained with 
GEOS-3 and SEASAT in non-repeat, or "general" orbits. This situation is 
likely to be reversed by future missions, where most satellites will be in 

repeating orbits. However, after the fuel of the stabilizing engines is 

exhausted, the satellites will gradually drift towards more "general" 
orbits under various disturbing influences (or may be maneuvered towards 
them while there is still some fuel left) yielding more suitable "master 
arcs" to rectify the mean sea surface, as long as their altimeters 

continue to operate for a substantial part of one apsidal period ( 2 u/ u) . 
In such orbits, most of the radial error of gravitational origin is 

observable in the crossover differences over that interval. 

According to (47) and (52), Ar( t j )- Ar( t^) in (92) is a function of 

the AC , so these may be estimated from crossover differences (Shum, 
nm a 

1982). The shorter period effects of the zonal errors are unobservable 

when the ground track repeats, and one would expect that also those 

produced by low order AC might be hard to observe (Wagner, ib. , 

ntnot 

1984). The situation is somewhat better in non-repeating orbits, where 

all AC are estimable, in principle. In any case, the use of crossover 
nma 

differences can be an interesting way of employing altimetry to improve or 
calibrate existing gravity field models, in combination with ordinary 
satelite tracking data, terrestrial gravity measurements, etc. 


6.3 Direct Mapping of Ocean Tides 


Because of the convenient structure of the orbit error, "frozen", 
repeat orbits provide the best data for studying tides, particularly in 
the form of overlap differences. According to expressions (85), ( 86 ), 

(87) and ( 88 ), the tide-related signal is the sum of At and Ar^, both 
of which, in the case of a repeating ground track, have exactly the same 
temporal frequencies. The overlap differences can be modelled using the 
observation equation 


*M*i» e J?) = ^ £ A? igT B i (d»,A)){-A (t y ) cos(a) T(central) t -8 7 ) 


y=l i Bx( central 


t(band) 


- B(ty) sin( t>J T(central) t y-B f) + A k ^cosMty + B^inMty 
t( band) 


+ C k tyCosMty + D k tySiriMty + E k tyCosMty 


+ F k t?,sinMt Y + i r t (t Y -t k )^ } (-1)^ Y_1 ^ , 
k Y ’ T j=0 ' T Y 


(93) 


where y=l corresponds to tj and y=2 to t . The terms of the sum "2" 

(X) ■ - ^ T ■ . . 

represent the "lumped tide" At^+Ar , with ^ T ^j }an( jj and ® T (b an d) 2 * ven 

by (70). In what follows, regardless of ty, the instant t^ always belongs 

to the first repeat of the ground track, thus A k "Aj , B k ^=Bj , etc. So all 

the overlap differences are with respect to the first arc , which should 

not be shorter than one repeat period T rc . The base functions B^ for 
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parameterizing the "lumped tide" in space can be anything convenient, 
including spherical harmonics if the estimate is world-wide. The velocity 
potential functions (see expression (78) et seq.) are orthogonal and 
complete over the ocean basins, so they can represent the "lumped tide". 


Solving a redundant system of observation equations based on (93) by 

the least squares method, one can separate the non-periodical part 

of Ar from the "lumped tide". Then, this "lumped tide" may be corrected 

iteratively (for example) to obtain A? q . To start with, one can guess 
(T) 

that Ar + A?^ is identical to A? q , making in this way a first 

estimate of the yielding of the solid Earth 6" and the perturbation of 

U? o ) 

the potential V . This requires a numerical spherical harmonic 

analysis^ ^ of that guess of A? and then use of the resulting At 

o nmaBi 

according to (73) and (86) to get the corresponding values of 6' and 

( T) 

Ar . Correcting the "lumped tide" by substracting these values will 

result in a new estimate of At , which can be corrected likewise in turn, 

(T) 

etc. Because Ar and S' are probably small compared to At Q , the 


d^It is possible to do this numerical analysis, to very high degree and 
order, using fast algorithms like those described (including program 
^Listings) in (Colombo, 1981). Recently Tscherning et al. (1983) have 
compared the efficiency of these and other algorithms, all of which 
approximate integration over the sphere with numerical quadratures on 
grided data. For sufficiently fine grids, computing the coefficients in 
this way, reconstituting the data from the computed coefficients, finding 
the differences between the reconstituted and the original data, and 
iterating this procedure, can be shown to be equivalent, in the limit, to 
a full least squares adjustment of the coefficients. The first iteration 
corresponds to the well-known practice of "pre-whitening" in ordinary 
Fourier analysis, and may increase the number of significant figures 
recovered considerably. To avoid numerical problems, the coefficients 

calculated in these large analyses are always fully normalized. 


iterations are likely to converge. As <0.1, 6' is one order of 

(At ) 

magnitude less than 5 q , and so is the radial effect of AV on the 

adjusted orbit, or epheraeris (5 cm at most for STARLETTE (1000 ra altitude) 
according to Marsh et al. ib. (1983)). The terms in Ar v ' can be 

separated in two groups: those with periods of the order of months, which 
are filtered out together with the non-tidal orbit error, and those with 
periods of less than T rc , that remain in the "lumped tide". Of the 

latter, only terms whose frequencies are less than one cycle per day may 
have to be corrected for, as their amplitudes decrease quickly with the 

period. The process described here assumes that the Love numbers are 
adequately known, the same as the solid Earth tide 6, Moreover, the 

contribution to 5' and Ar^ of the tides in the polar regions not 

spanned by the ground track would have to be estimated relying on pre- 
existing tidal charts. The use of the 4^ for parameterizing the 

successive approximations to AC q , starting with A^, may be advantageous 
because of their "tide-like" nature, and their orthogonality over the 

oceans . 

Once Ar and 6" have been separated from ^ in a global way, a method 
like that of Cartwright and Alcock (1981) can be used for the regional 
refinement of the tidal map. Estes (ib. , 1980) has simulated a global 
estimation of tides from altimetry (with a somewhat simplistic error 
model) and his work suggests that a faithful picture may be obtained from 
a few months' worth of data. A satellite with a coverage as dense as the 
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one projected for POSEIDON (Figure 2) should help greatly to obtain very 
detailed charts. 

If the geocentric tide can be mapped reliably with altimetry, it 
may be possible to use the measurements of taken with deep-sea 

pressure gauges (Zetler,, 1980), to estimate 6" + 6 directly at the 
locations of those gauges. This may be helpful in studying ocean loading 
effects, the body tide, and the rheology of the Earth's interior. 

A quite conservative s ufficient condition for separating the 

(NTI (TV 

non-tidal orbital error Ar' ' from the tidal signal A?_+ Ar is that 

G 

their respective spectra should have no frequency in common. In a 

(NT) 

"frozen", repeat orbit the portion of Ar present in the overlap 

differences consists of the "quadratic once per revolution" oscillation 

and the polynomial terms in (52). Its spectrum is confined to two bands: 

the first extends, approximately, from 0 to 2 oj and the second from 

M-2w to 1H-2S) (u> ~ one cycle every 30 years for SEASAT). The tidal 

signal, as seen from the satellite according to (86), (88) and (89) has 

all its components at frequencies that differ from the harmonics of ui^by 

± to , ,v. As 0 and M are harmonics of u> , the sufficient 

x(strobed) re’ 

condition requires that 

2 1 “I < I “ T ( strobed) I (94) 
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so all the spectral lines of Ar( T ) fall outside the two main bands 
of Ar( NT ). However, this guarantees good separability only if the 
duration T dr of the data stream analyzed satisfies 

Tdr > I 2»/C I "-Ttstrobed) I -2 l^l ) I * < 95 > 

For all past altimeter missions, and for any contemplated future ones , 
Tj < 4 years, so the tidal components that can be resolved from the 
measurements of. a single satellite must have "strobed" periods not much 
longer than 4 years, assuming that the data were free from errors. As the 
period of 2u> must be of the order of one decade, all those "resolvable" 
components should satisfy (95) and be, therefore, separable from the 

resonant terms in (93). 

In the case of ordinary non-repeating orbits where the perigee 
precesses once or more per year, expression (90) shows that the orbit 
error must resemble a periodical up and down motion with its main 
components at frequencies 0, u> and 2u>, and with spatial wavelengths as 

small as 1000 km. These components can be seen as additional "tides", 

which would have to be included in the model, greatly increasing the 
number of parameters to be adjusted; otherwise, they may bias the 

estimates of the real tides. 

All attempts at mapping with GEOS-3 data, originally thought to 

be a reasonably straightforward operation (Zetler and Maul, 1971), have 
turned out rather disappointing results. The poor distribution in time 
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and in space of these measurements (collected over disjoint intervals of 
no more than 20 minutes, because of the lack of an onboard data storage 
device), the somewhat low accuracy of the instrument and the non-repeating 
nature of the orbit, all have conspired to frustrate those efforts (Maul 
and Yanaway (1978), Bretreger (1979), Won and Miller (1979), Masters et 
alt (1979) and Coleman (ib., 1981)). The first estimates of the M2 tide 
showing clear agreement with the amplitudes and phases observed at 
maritime stations or plotted in existing tidal maps have been those of 
Cartwright and Alcock (ib. , 1981) for the northeastern Atlantic, of 

Mazzega (1983) for the Indian Ocean, and of MazZega (1984, ib) for the 
whole world, all based on measurements from the last month of SEASAT. 
Although these workers have taken no advantage of the structure of the 
orbit error (treated by them like white noise), their results are better 
probably because of the special nature of Ar along the repeating ground 
track, while the accuracy and the spatial coverage of the data must have 
helped also, even when the total observing period was short. 


7.0 CONCLUSIONS 


The nature of the orbit determines that of the radial orbit error. 
There is a geographical component in this error, due to the gravitational 
field model, that is largest in the case of "frozen”, repeat orbits. 
Because of the precise repetition of the ground track every few days, this 
is the kind of orbit likely to be used in most future missions, and has 
been used already for SEASAT. 

In the differences in altimetric heights along the colinear passes of 
a repeating ground track, the radial error left consists of a slowly 
increasing oscillation plus a low degree polynomial (if the arc is about a 
week long). This error, produced mostly by resonances of zonal origin, is 
easy to model and to separate from oceanic variations because it has a 
much longer wavelength. Its main contribution to the power spectrum of 
the total error should be confined to very low frequencies and to a narrow 
peak at about one cycle per revolution; this seems to agree well with 
experience (Marsh and Williamson, ib. , 1980) . It is important to have 

accurate and dense tracking and also good models of the non-gravitational 
forces. The use of incorrect models of these forces may increase 
significantly the resonant part of the radial error. It is easy to deal 
with the long period effects of the errors in non-gravitational models, 
because they belong mostly to the "resonant" part of Ar. The shorter 
period effects of these errors, on the other hand, cannot be treated 
adequately with the type of analytical theory used here. Therefore, the 

v: l 52 


present conclusions (insofar as they are based on that theory) are 
reliable only if the models are good enough to make such effects quite 


small. 


To separate the main tidal components from each other and from the 
mean sea surface, heliosynchronous and high inclination orbits should be 
avoided. However, for studying other time-varying features, using repeat 
heliosynchronous orbits could make the effects of non-rgravitational force 
model errors more tractable, particularly those caused by the complex and 
often changing shapes of the spacecraft. If the attitudes of the various 
parts of the satellite (particularly the solar panels on which most of 
those forces act) are stabilized with respect to the Sun by making their 
orientations in inertial space functions only of the spacecraft's position 
relative to the plane of the ecliptic, then the pressures of solar 
radiation and of Earth's re-radiation, and the errors in their models, 
will be functions of that position as well (except for small variations in 
the Earth's albedo, the slow annual changes in solar declination and in 
apparent luminosity of the Sun with distance as the Earth moves in its 
elliptical orbit, etc.) Under the assumptions, over long periods of time, 

i the errors in the computed values of these forces will repeat themselves 

■I. : ■ 7 _ ^ • 

} almost exactly once per revolution. Their contributions to the radial 

I error will consist of a resonant part (an extreme example of which seems 

F ■ , : 7 : .... 7 : : ■ : ' : 

jjj . ■ ■ ;• ; . 

to be present in one of the graphs in Figure 4), and of a part that 

| repeats along the ground track and thus largely disappears from overlap 

f differences. As the atmospheric helium bulge follows the Sun, similar 



considerations would apply to drag model errors. Departures from this 
well-behaved pattern due to changes in solar radiation Intensity and 
magnetic flux could be reduced by such a simple procedure as the daily 
adjustment of the drag and reflectivity coefficients (reality, of course, 
can be wonderfully messy, so these ideas are to be taken with a grain of 
salt). 

In a "frozen", repeat orbit, part of the radial error is "lumped 
together", in the altimeter measurements, with the error in the geocentric 
tide, but they can be separated using Love numbers. 

"Frozen", repeat orbits are best for studying temporal variations 
(like tides) by analyzing collnear pass differences, because the radial 
error left in them is easy to filter out; however, since a significant 
part of Ar is also unobservable in crossover differences, these orbits are 
worst for modelling the mean sea surface by the "bias + tilt" method. 

In the crossover differences of ordinary, non-repeating orbits, where 
the mean perigee is not "frozen", the radial error can be observed and 

corrected better than in the repeating case, by the "bias + tilt" method, 
provided the times of the ascending and descending passes at the 

crossovers are distributed over a substantial part of one apsidal cycle. 
In this sense, non-repeat orbits are best for mapping the mean sea 

surface. The same applies to the use of crossover differences for 

calibrating or correcting existing gravity field models. 
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Improving the force models, the gravitational one in particular, is 
essential to the determination of the stationary sea surface topography 
w Q , which requires both precise knowledge of the orbit and of the geoid. 

The difference between the geocentric tide (observed with satellite- 
borne altimeters) and the "surface-bottom” tide (sensed with deep-sea 
pressure gauges) is the sum of the body tide 6 and the crustal 
bending Once long series of accurate measurements from altimeters and 

pressure gauges are available, it may become possible to map these 
vertical movements directly, across the oceans. 
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